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iiBSTRACT 

In this thesis a two lane bridge made of pre stressed 
concrete has been optimized for minimum cost. The cost of 
bridge deck consists of the cost of longitudinal and 
transverse beams. The analysis of the deck has been via 
equivalent orthotropic plate and individual meriibers forming 
the deck have been designed by working load method for loa.ds 
and stresses specified by Indian Roads Congress. 

In this thesis two commonly used sections in the bridge 
decks have been taken up. The sections . are T-se ction and 
Hollow Box section. The non-linear programming problem of 
finding the minim'um cost design has been solved by sequential 
unconstrained minimization technique using interior penalty 
function. Davidon-Fletcher-Powell algorithm has been used 
for unconstrained minimization, while Fibonacci method is 
employed for linear search. 

Graphs has been generated to help in calcula.ting the 
maximum transverse moment in bridge deck for few classes of 
loadings of Indian Roads Congress. A parametric study has 
been performed with number of longitudinal beams in a deck 
changing from 5 to 9. The spans considered are 24m and 3om. 
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INTRODUCTION 

1. 1 Introduction 

Simply supported pre stressed concrete "bridges for 
highways consist of a momber of longitudinal beajii.s which are 
hinged to the support at one end and resting on rollers at 
the other. These longitudinal beams are made monolithic with 
the transverse beams which are also of pre stressed concrete. 

The nximber of transverse beams are usually predetermined. With 
such a configuration of beams in a bridge deck, it can be 
idealized as a gridwork. The analysis of gridwork subject to 
transverse eccentric loads by transforming it to an equivalent 
orthotropic plate was pioneered by Guyon [l]. It was further 
developed by Massonnet [2]. The formal technique of analysis 
and design of such a gri'dwork via equivalent orthotropic plate 
was fully documented by Rowe [ 3 ]. Sven today it is a classical 
approach for Bridge Designers. 

For a developing coirntry like India where the transport- 
ation system is yet to be perfected, highways and hence 
bridges will assune great importance. Prestressing techniques 
have been understood well by now. m number of small factories 
exist throughout the country where pre-cast sections can be 
made and transported in convenient si?;es to the site and 
then post-tensioned in situ^. In the present work few commonly 
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used sections for pre-cast and post-tensioned concrete bridge 
decks have been optimized for minimum cost. The loads acting^ 
on the bridge are as per the criteria laid down by Indian 
Roads Congress (IRG)[5]. The working load method has been 
adopted for design. Optimum design is arrived at by sequentia 
unconstrained minimization technique [6] using interior 
penalty function. The David on-Pletcher-Powell algorithm has 
been used for unconstrained minimization and Pibonacol method 
is employed for linear search. 

1.2 Siunmarv of Previous Work 

There has been sufficient published work on the analysis 
of bridge decks [3,4,7J. However, little work on optimum 
design of bridge decks is available. A completely automated 
design of a three span bridge of prestressed concrete has 
been developed by Sawko and ¥illooofc [8] to satisfy the 
lo3.(iiiig P6Q^ijiiP0in6n "ts « of Bpitisii 

The handbook [9] prepared by Structural Engineering 
Research Centre (SERC) gives optimum standard sections for 
the loading proposed by IRC. The sections considered are 
inverted T-beam, Hollow Box, composite I section and T 
section for the spans ranging from 12m to 36m at an interval 
1. 5m. The clear carriage way taken is 7.34m. The number of 
beams considered incase of T-section bridge is 5 , 6, 7 , 8 
and 9. The cost of these bridges is compared. A bridge 
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having five beams turn out to be of minimum cost. In the 
case of Hollow Box section bridge, the number of beams 
considered is 9 and 12. The cost of the Hollow Box bridge 
having nine beams turns out to be the minimum. While calcu- 
lating the cost of bridge deck the cost of prsstressing cables 
in transverse beams are not taken into account. Moreover 
the minimum cost is arrived at by a brute force method of 
comparing the cost of bridges with beams increa,sing in depth 
by one centimet€T~. The distribution coefficients for 
longitudinal moments corresponding to unit point load are 
read from graphs given in reference [3] for various standard 
reference points on the deck. 

The optimum design of T-beam bridges subject to 70R. 
loading of IRC was carried out by Raman et al [lo] satisfying 
the stress condition of IS-1343. The nonlinear optimum design 
problem has been solved by linearizing the functions through 
Taylor series. The method used for distribution of longi- 
tudinal moment is due to Oourbon [ll], which is a very 
approximate method. 

Kumarasamy [l2] minimized the total cost of a T-beam 
bridge having several simply supported spans using seq_uential 
linear programming technique. The cost function was minimized 
subject to constraints stipulated by IRC. The variables in 
the optimum design problem was number of beams, number of 
bays, nmber of prestressing cables besides the other 
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variables which define the geometry of the section. Here 
also the cost of transverse beams was completely neglected. 
The coefficients of distribution of longitudinal moment 
corresponding to unit point loads was supplied at discrete 
points. The graphs have been drawn to directly get the 
coefficients for distribution of longitudinal moments for 
two lane bridge for 70R, class AA and class A loe^ds. 

The T-beam. sections used in prestressed concrete bridge 
were optimized by Kapoor and Bandhopadhyay [13] for a simple 
loading. The working load method, ultimate load method and 
reliability based method were used to optimize the cost of 
the section. It was shown that reliability based method 
gives the least cost design. However, the investigation 
suffers with the drawback that practical loadings have not 
been considered. 

1 , 3 The Present Work 

A two lane bridge with foot path on both sides has been 
optimized to obtain the minimum cost of the entire deck 
system. The following two commonly used sections in bridge 
decks have been chosen for study in the present work. 

( i) T-section 
(ii) Hollow Box section 

The cost includes the cost of longitudinal beams, 
transverse beams and the prestressing cables. Equivalent 
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orthotropic plate analysis is used to determine the worst 
loaded beam in both longitudinal and transverse directions. 

All the longitudinal beams are of same cross-sectional area, 
shape and size. The transverse beams are also identical 
in ' detail. 

A brief description of orthotropic pla,te analysis of 
bridge decks is given in Chapter II. The gra^phs developed 
for coefficient of distribution of transverse moments are 
also presented in this chapter. The graphs correspond to 
70R and class AA loading placed to cause maximum moment 
in the transverse beams. 

The working load method of design is described in brief 
in Chapter III. The formulation of the optimum design 
problem and the solution technique is also described in 
short in this chapter. 

In the last chapter the results of two problems, one 
for each type of section, are presented. It is followed by 
discussion of the results. 

It is concluded that T-section bridge deck is economical 
in comparison to Hollow Box section bridge docks for the 
spans studied. The parametric study on these bridge decks 
with increasing number of longitudinal beams reveal tha.t 
the cost of bridge deck increases with increase of 
longitudinal beams. The use of 7 number of longitudinal beams 
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is recommended as the slight increase in cost of deck is 
expected to he w el- compensated in handling the reduced weight 
of the modules. 



GH/J>TER II 


ANALYSIS OF "BRinGE DECKS 
^ 1 Introduction 

There are various methods available for analysis of 
bridge decks [5,4,7,ll]. The method adopted in the present 
work is the one which converts the gridwork of bridge decks 
to an equivalent orthotropic plate. This approach seems to 
be commonly used by designers. It uses a single set of 
distribution coefficients for two extreme cases of a 
no-torsion grillage and a full torsion slab to find the 
distribution of behavior pararceter e.g. forces and deflections 
etc at any point of any type of the bridge. The governing 
differential equation of flexure for orthotropic plate and 
the boundary conditions, the equivalent orthotropic plate 
has to satisfy, is given in the following section. The 
different parameters involved in the solution of the 
governing differential equations are discussed alongwith 
definitions of coefficient for distribution of longitudinal 
and transverse moments. The procedure of finding the design 
moment for the longitudina,! beama as well as the transverse 
beams are discussed in the subsequent sections. The chapter 
ends with the section describing the graphs generated in the 
present work to help in calculating the transverse moments. 
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2 , 2 G-overnirLie: Equations of an Orthotroipic P late 

The theory of orthotropic plate is described in detail 
in reference 14. For completeness of presentation the 
equations in brief are given; 

The stress-strain relationship is given by 

“x = 


I I* 

a = E e + E s 

y y y x 


2.1 


T = Gv 

xy *xy 


Where e'" = = 




^ ' ^21 2 , 2 

'x y *^y 


assuming plane sections remain plane after bending^ the 
strains in terms of deflection ■w(x,y) are given by 


0 _ 2 




£ = -Z 

y 


3y 


2 * 


Y 


xy 


= -2z 


ax3y 


2.3 


The corresponding stress components are 

o = -z(E -5-7, + E -jr—,) 

X «x2 av^ 


X 


a = -z(B„ + E 


'y 3 y" 


ay‘ 

“ax' 


^ = .2Gz 

xy 6xay 


2.4 


The bending and twisting moments are given by 
h 

3^w -f. a^w 

1 3y 

2 


M = 

X 


XI 

hi 


dz = -(Dx g ^ 
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h 


M 


y = f ^ a e 4z = -(By + B, 

J- -2 


k 
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and = 

J . 

where Dx 

" 'Z 

= X 

__ 




a.:T' 




1 Zx‘ 


... 2.5 


1 z dz = 2Dxy 
xy '’2) xay 


E h3 

; Dy = -t— 


12 ' " 12 ’ 


Dzy = 


&h^ 

12 


. . 2.6 


By using moment equilibrium of the element, the relationship 


for shear forces Qx and 4y is given as 


^ M. 


■a M 


Qx = -r—- + ^ 

ax By 

a M a 

Qy ^ J 

ay ax 


_ _ _^(Dx + H— 

Tr a X^^^ -W 2 ^ -U ^ / 


a x‘ 


w 


= - ----(Dy ^-T + H- 
By ©y^ 


ay‘ 

a^x'‘ 


where H = D-, + 2D 


xy 


...2.7 


..* 2.-8 


By considering the static equilibrium of vertical forces 
acting on the element shown in Pig. 2,1, we get 

jw__ . T,— ,,. 2,9 


*^4* 4" 

a W , OTI ® ' 

Dx T + 2H 


By = 4 


ai? ' “ax=ay= ^ 3y4 

The equation (2.9) is the governing differential 
equation for flexure of ortho tropic plate where n is a 
measure of torsional rigidity of the plate. 


2 . 3 Bridge Deck and Squivalent Orthotronic Plate 

A bridge deck has equally spaced beams in both longitudinal 
and transverse directions. In general, the stiffness and 
spacin^s of these beams are different in different directions. 
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Thus if the bridge deck has to be transformed to an ortho- 
tropic plate, the elastic parameters relating to the 
substitute system are to be assumed as continuously 
distributed and their respective values are to be defined per 
unit length. 

Consider a grillage consisting of two systems of equally 
spaced parallel beams and loaded perpendicularly to x-y plane 
as shown in Figure 2.2. If I and J are flexural moments of 
inertia of individual beams along x-axis (longitudinal 
direction) and y-axis (transverse direction) respectively. 
Then 


Dx 




... 2.10 


where B is young's modulus, a^^ and b^ are spacing of beams 
in X and y directions respectively. Thus Dx and Dy represent 
flexural rigidity of equivalent system for bending in x and y 
directions respectively. 

Again if I and J are torsional moments of inertia 
of individual beams along x-axis and y-axis respectively, 


GI 

Y ^ ® 

^x - ^ ’ 


^y " 


. 2.11 


Y and y torsional rigidity in x and y directions 

X y 

respectively. 

The bending and twisting moments can now be written as 


^x = 


- Dx- 


x^ 


M = 

y 


- Dy 


3^w 


ay 


2 r 


M 


xy 


= Y 


and H 


3^w 
X a x'ay 


yx “"^y axay 


... 2.12 
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Considering the vertical equilibrium of an element, the 


governing differential equation of flexure comes out as 


or Dx 


B x4 


“b x4 


. / . X T, 3^w 


, OTT w o w 


..2.13 


. .2.14 


where 2H = y + Y 
X y 


. . .2.15 


and it represents the measure of torsional rigidity of the 
equivalent orthotropic plate. 

The equations (2,9) and (2.14) are of sajne mathematical 
form. If Poisson's ratio (y = 0.15 for concrete) is 
neglected the expressions for M^, M are same for equivalent 

y 

orthotropic plate and the true orthotropic plate. If Qy is 


defined as 


Qy = - DyS 


1 / \ 3 ^w 

+ Yy) ^x^3y 


. . . 2.16 


the expression for M and M are also identica.1 for the 

xy yx 

two systems and the value of approximately becomes zero. 
By comparing equations (2.15) and (2.8) 


+ Yy = 2H = 4Dxy 


. . .2.17 


2. 4 Boimdary Conditions 

The bridge is simply supported along the two opposite 
edges at X = 0 and x = B, and is free at the remaining two 
edges y = +b. 
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At the supported edges the deflection is obviously zero 
for non- sinking supports i.e, 

w = 0 at X = 0 and L 


md the bending moment M =0, i.e. 


...2.18 




^X‘ 


0 at X = 0 and L 


Along the free edges bending moments M is zero, i.e. 


y 


9 

■ “ ' i H 

sy 


= 0 at y = +b 


and the shear is also zero i.e. 

Qy = Dy + 2Dxy = 0 ; at y = +b . . .2.19 


2.5 Dimensionless Parameters 

In all practical cases the value of torsional rigidity 

given by the sima (Yv + Y ) enclosed between two limits 

X y 

corresponding to the simple grid of beams on the one hand 
and the isotropic slab on the other. 

Putting 


Y^ + Yy = 2a V'DxDy 


. 2.20 


the equation of elastic surface (in terms of deflection w) 
in the simplified form is 
w 


S'* 


^ 4- 

© -W 


•.^4 

d w 




2.21 


and for the shear force Qy 

Qy = - By - cc\rnSi gl— ...2.22 
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In the last two equations the dimensionless parameter a 
is characteristic for the resistance in torsion. The range 
of variation of a is limited hy the values 0 and 1. For a = 0 
the system reverts back to a very simple grid of beams with 
no torsional rigidity while for isotropic plate a = 1. From 
Equation (2.20) we have 


a = 


Y + Y 

y 

2Vl)xi)y 


. . .2.25 


Equation (2.25) can also be written as follows: 


a 


Y + Y 

Is ^ 

2fDxDy 


ZS'iTT 


... 2. 24 


where 



= Torsional inertia per unit width of bridge 
= Torsional inertia per \mit length of bridge 
= Moment of inertia per unit width of bridge 
= Moment of inertia per unit length of bridge 


The determination of and is a. tedious job. For 
a T-section which consists of three flat plate elements, 
Figure 2.3, or is given in reference 5 

I or J = k,(2aj^*2b. + k-(2a„)^-2b„ + ik_(2a,)^2b„ 

0 o ii i d d d 3 j j 


...2.25 
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where suffixes refer to rectangle 1, 2 and 5. The value of 
b corresponds to the longer side of the rectangles. The 
torsional constants k used in Equation (2.25) are given by [l5] 



M ^ 

'^i m=l,3,5. 



tanh 



. . .2.26 



The width of the equivalent orthotropic plate tnms 
out to be equal to the width of the bridge, because all the 
longitudinal beams considered in the present work are 
identical which generally happens to be the case in all 
practical designs. 

a and 9 given by equa.tions (2,24) and (2.29) respectively 



'X 
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i 2 W 
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fig. 2»5 s Bcceatrio staasaidal load m 
bridge deck 
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are for homogeneous cross-sections. In case of prestressed 
concrete sections where no tensile stresses are allowed, the 
cross-sections can he considered to be homogeneous. IRC 
does not permit any tensile stresses to develop in these 
sections. Hence the sections considered in the present work 
are homogeneous in this context. 


2 . 6 Principal Coefficient of Distribution of Longitudinal 
Moments 

The load acting on the equivalent orthotropic plate 
is ass-umed to have a sinusoidal va.riation along the span of 
bridge, or 

p(x) = sin ^ ...2.30 

and the load is assumed to be ^^cting at an eccentricity 
’e' (Fig. 2.5) in the y-direction. Due to this load the 
deflection can be assumed as 


n X 


..2.31 


w(x,y) = W(y) sin 

If the load defined by eqmtion (2.30) is distributed 
over the entire width 2b of the system, the intensity of 
this load, uniform in the y-direction becomes 

p (x) = p sin where = — 


o'" -"0 1 

and deflection is of the form 


2b 


w (x) = ¥ sin 
0 ' o 


2£X 

T, 


2.32 
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The ratio of the two deflections, as produced by the 
line-load and hy the distributed load p^ respectively 
is called the prinaipal coefficient of distribution of 
longitudinal moments, i.e. 


o o 

The value of K(y) will depend on the parameters 9 and a, 
relative eccentricity e/b, and ratio y/b the rele^tive 
ordinate of the section considered. Moreover, K(y) gives the 
deflection at a point y/b for a unit sinusoidal load acting 
B.t an eccentricity e/b. 

Massonnet [l6] has solved the equation 2.14 for a. load 

The system analysed is a grid with no 


p(x) = sin ^ 


nx 


torsional resistance. The deflected surface is given as 

^1^ X 1 / 

" = "51^ 1 Sinh?2>,V- -sin^S^bV 2ooshA(b+y) ooSA(b+y) 


[sinh 2?vb cos?k(b+e) cosh?\(b-e) - sin2Pib coshA(b+e)' 


cosJ\(b-e)] + [coshA(b+y) sin;v(b+y) + sinhA(b+y) cos?v(b+y)] 
^sinh2Ab [sin(b+e) cosh;s(b-e) - cosA(b+e) sinhA(b-e)3 ; 

+ sin2Ab[sinh?>(b+e) cosA(b-e) - coshA(b+e) sinA(b-e)'|^ 

...2.34 

’Dx 


where A = 4 


71 

in 


Dy 


tiQ 

bV2 


.2.35 


If the load P]_sin jj- is loniformly distributed over the 
entire width 2b the system will deflect to a cylinder as 
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given ly equation 2.37. By dividing equation 2*39 "by 2.37 
the principal coefficient of lateral distribution is obtained 
as 

*^0 = slnh^2^b~- 3ln^2At ( cosA(b+y) [sinh2Ab - 

cosS\(b+e) coshA(b-e) - sin27*b coshT'Cb+e) cos?\(b-e)] 

+ [coshA(b+y) sin^(b+y) + sin^?k(b+y) cosA(b+y)] 
j^sinh2?sb[sin>k( b+e ) ccshh(b-e) - cos>>(b+e) sinhA(b-e)J 
+ sin2^vb[sinhA(b+e ) cos/^(b-e) - cosh>\(b+e) sinX(b-e)]^| 

.. . 2.36 


The subscript ’O' denotes that the formula "holds true 
for a = 0. The equation (2.36) is valid for the portion of 
the width of the bridge where the ordinate y is less than 
the eccentricity, e, of the load. In case we are on a point 
in the bridge where y is greater than e, the sign of y and e 
in the above equation get changed. 

For a = 1 the corresponding coefficient can be 
calculated using the formula, of Guyon [l] given as 

Kf = [(ccosho + sinhcr)cosh ©"X- ©"X 

R Q Qu/ 

4. 2 — !l 4. ^1 2 37 

3sin,h .c cosha - a 3sinh0 coshc + 0 ... .--i 


where 



(acosha - sinho) cosh©q) - ©(psinhc sinh©q) 

(acoshc - sinha) cosh©y - ©ysinhosinhQw . 2. 38(a) 
(2sinh0 + acoshc ) sinh©9 - ©9Sinha cosh©9 
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Qy = (2sinha + crcoshcs’)sinh9y - 94»sinha ccsh©^ 

y = ^ ...2.58(b) 

ff=©Ti; 

lor a practical case the value of a lies between zero and one. 
So to find E an interpolation formula 

+ (K^ - EQ)f5 ...2.59 

is used. The value of given by equation (2.59) is increased 
by 10 percent to make the theore^tical results correspond to 
actual experimental results [5J. 

For a bridge of known dimensions Q and a can be easily 
calculated. E^ and can be found for a given value of 9 
and for different values of eccentricity of load, e/b ranging 
from -b to b at intervals of b/4, and for different discrete 
locations on the bridge, y/b ranging from b to b at intervals 
of b/4. These values of and are stored in two (5x9) 
matrices. Using Maxwell's theorem the two (5x9) matrices 
ca.n be expanded to (9 x 9) matrices to include the effect at 
other points of the bridge corresponding to y/b ranging from 
-b to o at interva,ls of b/4. The nine discrete points 
(“b, - ~ "i’ hence forth referred 

to as reference stations. 
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2 Maximm Longitudinal Moments 

Maximxim longitudinal moment can be calculated from the 

two m8.t rices [K^] and [K^] for any practical case, using the 

interpolation formula given by eq_uation (2.39). 

For the spans under considere.tion (24m and 30m) it ha.s 

been established [l2] that 70R loading of IRC gives maximum 

longitudinal moments. The maximm eccentricity of load is 

at a distance of 1.2 metre from the edge of the foot path. 

The actual load position for worst effect is shown in 

Fig. 2.6. As the wheel loads do not coincide with the 

reference stations, so these loads are repla.ced by equivalent 

loads at the reference stations. This is done to facilitate 

the use of matrices [K^] and which have been generated for 

these reference stations. The equivalent loads at the 

standard reference stations are the reactions of the simple 

beam of span b/ 4 [3]. The loa.ds at the reference stations 

4 

are thus >>41 where ^ = 4.0. ?\’s are 

known as weighting factors. The maximum elemorit of the 
vector obtained by 

|C?v]^[K^] + - K^] 2.40 

will indicate the location and the value of maximum 

coefficient of distribution of longitudinal mom.ent. Since a 

longitudinal beam may not lie under this reference station, 

this coefficient (K ) is interpolated for the nearest 

amax 
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beam centreline, a.ssuining a lineair variation. Thus 


M = I.IM K 

xmax xmean aina.x 


...2.41 


where M is avera.ge moment per beam of the bridge, equal 

to the applied moment divided by number of beams. M is 

the bending moment for which each longitudinal bean will be 
designed.. 


2 . 8 Coefficient for Distribution of Transverse Mor a ents 

The transverse moment M is given by equation (2.12). 

y 

If w as given by equation (2.35) is put in equation (2.12) 
we get M in the form 

y 


M = 


, ,T • mx X 

bHm sxn 


. . . 2.42 


where p „ is a coefficient dependent on 9, a, e/b and y/b, 
my 

a.nd Hm is the amplitude of the load for various terns in 
Fourier's series. In the two limiting cases p^ and p^ are 
defined for a = o and 1 respectively in the same way as 
and were defined. We have thus [l6] 


Tz 'n;0( sinh^2Ab - sin^2/?ib) 


^2sinh>k,( b+y ) sin?v(b+y) 


[sinh2Ab cos^(b+e) cosh?\(b-e) - sin2Ab cosh>s(b+e) • 
cosA(b-e)] + [ cosh?v(b+y) sinA(b+y ) - sinh?v( b+y )cos?v( b+y) ] 
sinh2Ab[sin?v(b+e) cosh?k(b-e) - cos^(b+e) sinh>\(b-e)] 

+ sin2Ab [sinhX( b+e ) cosA(b-e) - cosh?v(b4-e) sin>\(b-e) i 

...2.43 ' 
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and 


^1 ” 4 crsinh^cr 

{ I [( 1 - y)gcosh.g - ( l+y)sinh.g jcoshQy - Q¥(l-v)sinhg sinfaQ ^ii/} 

\ X 3 + V ) s inhg c o s hg - ( l-i/)g 

•|'[( l-i>)gcosh.g - ( 3 + v) 8 inhg ]cosh 0 cp - ( l-:w) sinhg- 09 - sinlaQ 9 | 

r ( l-^)gcoshg 4- 2sinhg ]slnh0 V - ( l-if) sinh g* 04* cdshG i^ii . 1 
( 3 +i^)sinha cosha + (l-v)g J 

(l-i>)gcoshg sinh 09 - ( l-v)sinh.g* 09 - cosh 09 ^ 

+ [( 1 - J/)gcoshg - ( l+l,>)sinb.g ]cosh0“5<. - ( l-a/)s inhg- sin h©"Xj 


...2.44 

where u is poisson's ratio. 

Again for a practical ca,se a lies hetween zero and ione. 
So is interpolated according to the formula 

— Pq ( p-^ P|^ )V^ ... 2,45 

2,9 Maximum Transverse Moments 

The wheel position for 70R load to cause maximum 
transverse "bending moment at the midspan is shown in fig. 2,7. 
Rowe has shown [ 3 ] tha.t the maximum transverse moment occurs 
at X = and y = 0. In the present work transverse "beams 
are provided at supports, quarter points and midspan. All 
transverse "beams are designed and detailed identical to that 
at midspan. Due to these reasons, Rowe gives graphs for 
and at reference station 0. 
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Coming back to the equation (2.42), is the value of 

as given by 2,45 for flexural parameter m9. However, 
since only the centre of the bridge need be considered, the 
appropriate values of x in the equation (2.42) is x = 

niTC X 

hence sin —jj- is alternately plus and minus one for the odd 
terms of the series and zero for even terms of the series. 

The value of in this equation can be determined from 

the influence curves for and as given by equations 
(2.43) and (2.44) respectively for flexural parameter ©, 3©, 

5© etc. summation of the 

ordinates in the influence curves for and p^ respectively, 
the ordinates being measured under the specified position 
of wheels in the transverse section, 

Hm is the amplitude of m"^^ term in the Fourier series 
for the load. The amplitude of the Fourier series for TOR 
loading as given in [4] is 


Hm = 




mitu. 


mmu. 


mmu 


+ sin — : 


+ sin 


7 


} 


.. .2.54 


L I ■ • * • ’ ----- 

where u’s are distances of the respective axles from one 
support. Thus the maximum transverse moment can be written 
as 


M 


ymax 


21 

L 


r ^ 7 Tiu 7 

1 ^ ^ sin -|^) - 

7 


3tcu 


511 u. 


1^5©^ Si ^n *L 


sin 


h' 


1 


n “ 


...2.55 


n 


) 
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where 2. ...2,56 

^^ymax design moment for transverse heams. 

2.10 Graph for 2 M-q and 2 

It is obvious from equation (2.55) and (2.56) that 
computing the design transverse moment, as a tedious 

job. The determination of 2. and 2! involves drawing 

six influence curves, three each for p^ and p^, for 9, 39 
and 59 from graphs given by Eowe [3]. The ordinates at 
specified load position is marked on curves and summed up 
for a value of 9 to get 2 p^ and ^ p^^* Here an attempt 
has been made to draw the curve of p and "Zp, versus 9 
for TOR and class AA loading. The different graphs are as 
follows. 

(i) Figure 2.8 is for Tracked class AA load travelling with 
one train of load coinciding with longitudinal centre 
line of the bridge. 9 is from 0 to 5.0, 

(ii) Figure 2*9 is for Tracked class AA load travelling the 
bridge in such a way that the two train of loads are 
symmetrically placed about the longitudinal centre-line. 
After 9>2.25_ ^ !^q becomes negative. 

(iii) Fig, : 2. 10 is for TOR wheeled loads when one of the 

inner train of axle loads coincides with the longitudinal 
centre-line of the bridge. 
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(iv) 5lg. 2., 11 is for TOR wheeled loads when the train of 
axle loads are symrae trie ally placed about the longi- 
tudinal centre-line of bridge. Here also after 0 4.0, 

both "2. |i^ and '2 are negative. 

(v) Fig. 2.12 is same as Fig. 2„10* The difference is that 
in equation (2.44) vPoisson's ratio has been made zero, 

(vi) Fig* 2.13 is same as Fig, 2,ll except that V = 0 in 
this case. 

IBM 1800 was used to generate 250 points on these curves. 

The curves were drawn using IBM 1627 plotter. 



CHAPTER III 


DESIGR OP BRDIGB DECK 

3 . 1 Introduction 

The working load method has been used in designing the 
individual beams forming the deck. In fully prestressed beams 
cracks are not allowed to be developed in the section and 
material thus behaves as homogeneous and linearly elastic. 
Since the I.R. C. recommendation [17] does not allovj any 
tensile stress to be developed in the beam, the use of working 
load method for designing the beams is justified. 

The next section begins with a brief discussion on 
working load design. This is followed by statement of the 
problems taken up in the present work and the formulation of 
the optimum design problem. Chapter ends up with discussion 
on algorithms used for unconstrained minimization and linear 
search. 


5 ♦ 2 Working load Method of Design of Beans 

Let Pi be the prestressing force required at an 
eccentricity e for a prestressed concrete beam of cross- 
sectional area P., The section modulus for top fibre of beam 
be Z. , the section modulus for bottom fibre of the beam be 
Z-^ and applied livie load bending moment be while is 
dead load bending moment. 

The stress in bottom fibre at transfer condition is given 


by 


^i® ^i 

~z X~ 

^b ^ 


M 


bi 


...3.1 
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where is allowable stress in bottom fibre at transfer. 

The stress in top fibre at transfer condition is given by 

F.. M 


F e 

i i ■'''d 


ti 


. . .3.2 


where f , . is allowable stress in top fibre at transfer. 

"C 1 

Similarly the top fibre stress and bottom fibre stress at 

working condition are given by 

M, M 


- 4. jr n 2: -1. ^ 

■ + A 


and 


n^iS pFj_ 

-z"^-T 


'-'a 


. ..3.3 


. ..3.4 


"b "'b 

respectively. Where f^ is allowable stress in top fibre, 
f-j^ is allowable stress in bottom fibre at working conditions 
and p is loss in prestressing force. 

For a section of known concrete grade and sectional 
properties only F^^ and e are unknowns in the equations (3.1) 
to (3.4). But before F^^ and e are foimd, expressions for 
minimum section modulus required at top and bottom fibres in 
a beam will be developed. 

By multiplying equation (3.1) by p and substracting 


equation (3.4) from it 

Zy '^'^^bi ^b 

if equation (3.5) is taken as an equality the minimum section 


M M 

(1-p)^- + ■“ - A 


3.5 


modulus required at bottom fibre Z, is given by 


M -i- 

a 


bmin 


'i*bi 


. ..3.6 



Similarly from equation (3»2) and (3.3) the minimum 


section modulus required at top is given ly 


^ + ( 1 -T 1 )M^ 

Hmin f^ - 


. ...3.7 


The values of minimum eccentricity and corresponding 
prestressing force is worked out by satisfying equations 
(3.1) and (3.3), which gives 


f.Z, + pf, .-Z, - M - (l-p)M, 

.r, t t ' bi b a ' d. A 

A 


z z + z 

p _ — i ^ ^ ( — _ 

min A A ^f+.Z. + pf, .Z 


“a ^ ""d l_!t!t 

Ma-3-TI=;nM; 


. ..3.8 

) 


-t“t ' ''^bi^b 

...3.9 

Similarly the maximum eccentricity and the corresponding 
prestressing force is worked out by satisfying equations 


(3.2) and (3.4), which gives 


^t ^b + Z^ 


max 


T~ 


A 


“d - ^ti^t 


+ Xl-p )M^ + f^Z^~’+ pf+^Z 


■)p 


a 


■ti"^t 


“a + + ’'Pih 


and = 


P 


(Zh- + Zj 


-A 


. . .3.10 

. . .3.11 


The actual eccentricity provided is generally one which 
gives a maximum of O.lh as cover. ¥here h is the overall 
depth of the beam. This eccentricity is accepted if it lies 

between e . and e as found by equations (3.9) and (3.10) 

ruiTi ihslx 

respectively. So once the actual eccentricity e is known 
the prestressing force can be found by taking equation (3.4) 



... 3.12 


as an equality. Thus 


P. = 
1 


a d b h 

^(e + ^) 


The cables are laid for in a parabolic profile.' 

With P^, e and other variables defining A, Z_^, kno'wn 
checks are made to find the feasibility of the solution. The 
checks are 

(i) The section modulus provided for top fibre Z^ has to be 
more than given by equation (3.7). 

(ii) The section modulus provided for bottom fibre Z^ has 
to be more than given by equation (3.6). 

(iii) The four stress inequalities given by equations (3.1) to 
(3.4) have to be satisfied. 

(iv) The eccentricity provided has to be between the values 
obtained by equations (3.9) and (3.10). 

(v) The section has to be under-reinforced as given by 
IS-1343-1965 [19]. 

(vi) The ultimate applied bending moment given by 2.5M + 

1.5M^ has to be less than ultimate capacity of section. 


3 . 3 T ransverse Beam Desiisn. 

The dead load on the bridge deck can be taken as 
uniformly distributed on the entire span. Due to this 
symmetric load all the longitudinal beams will deflect by 
an equal amount and the transverse beams do not get loaded. 
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Thus only moment due to live load on deck will cause stresses 
to develop in the transverse beam. The transverse beams are 
designed as a member which is in unifom stress due to 
prestressing force. This is possible if the eccentricity 
is zero. The expression for the four stress inequalities 
become 



f 


bi 



Hit 


M 

^ f 
2tt ^ 


17 


M 


y 


"bt 


'b 


... 3.13 

. .. 3.14 

. . . 3.15 

... 3.16 


where F, is transverse prestressing force, A is area of 

X X 

transverse beam, Z. , and Z are section modulus for top and 

XX Du 

bottom fibre of transverse beam. The actual prestressing is 
given by 


M 'A, 
Y- t. 


" n . 


. . . 3.17 


The stress in bottom fibre at transfer condition as given 
by equation (3.13) and stress in top fibre 8„t working 
condition as given by equation (3.15) is checked, to ascertain 
if they are within permissible limits. The cables 
corresponding to F^ is laid straight. This is done so that 
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it does not cross the cables of longitudinal beajus. 


3. 4 Design for Under-Reinforced Section 

To pre'went failure of prestressed beams by crushing of 
concrete the clause 6.8.2.2(c) of IS-1363 [19] is effective. 


This is a’Chievted.:. by putting an upper limit to the steel 


to be used in the 


il 


t sw 


b’d Fc 


where /r = 

tsw ts 


section given by 
Fs 

— ^ 0.240 


.3.18 


A. = area of steel provided 

0 S 

^Hsf " Vi’s 

b = width of the flange 
b' = width of the web 

d = depth from the compression edge to the centre 
of the prestressing steel in the tensile zone, 
t = average thickness of flange of a flanged member 
For such beams the ultimate moment of resistance (M ) is 

U.S 

calculated as follows [l9] 


(i) if t d/4' 


A. • Fs 

Mus= W^^ts-^"-^tl-0.75^?|:jr^] 


. ..3.19 


(ii) if t d/4 


“us - O-'i'S 


A. -P 
t sw S I 

b”-d-Fc-' 


4 - O.TFc- t(b-b' )(d-0.5t) 


. . .3.20 








- .:?*>■■ ’f;-’ ^A 


\r’>V, ’'Jy^ri*"^ y^Xl ^SdS 


’■>'-:vrp:v?^4M|g 




C.G.0F SECTION 


C,G. OF PRESTRE^OiNG G 
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K = 1.0 lif tensioned steel is effectively bonded to concrete, 

liS 

and equal to 0,7, when steel is unbonded. 

Fs is ultimate tensile strength of prestressing steel, 

Fc is cube strength of concrete at 28 days. 

5 • 5 Optimum Design Problem 

3.5.1 Design Variables; As already mentioned earlier two 
types of cross-sections used in the longitudinal beams of 
bridge decks have been selected for study in the present work, , 
they are (i) T-section and (ii) Hollow Box section. 

(i) T-section 

Figure 3.1(a) shows the transverse cross-section of a 
bridge of 2- lane with foot path' on both sides. If number of 
beams used in the longitudinal direction is known (say N) 
and 2b is the total width of bridge, the width of top flange 
becomes automatically fixed given by WIF = ...3.21 

The Figure 3, 1(b) represents the idealized T-section 
used in the deck, the dimension which completely define the 
geometry of the section are XI to X5. The transverse beam 
is shown in Figure 3.1(c). The -prestressing force in 
longitudinal direction and transverse directions and their 
eccentricities can be found from these dimensions. 

The design vector for T-section bridge, therefore, is 
[X]^ = [XI, X2, X3, X4, X5] ...3.22 
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(ii) Hollow Box section. 

The Figure 3.2(a) shows the transverse section of a 
Bridge deck made up of number of Hollow Box section Beams, 

Figure 3.2(B) shows the idealized section of beans in 
longitudinal direction and Figure 3.3(c) is idealized section 
of Beams in transverse direction. The variable are XI to X4 
as shown. •0?op and Bottom width of flanges are fixed with 

the number of Beams predetermined. Here also the transverse 
and longitudinal prestressing forces and their eccentricities 
can Be found once the variables XI to X4 are assigned. 

Thus the design vector in this problem is 

[X]^ = [XI, X2, X3, X4] ...3.23 

3.5.2 Objective Function: The objective function (OBJ) for 
both problems are the cost of deck given By 

OBJ = Cost of concrete in longitudinal and transverse beams (Cl) 
+ Cost of pre stressing cables with its sheathing and 
anchorage in longitudinal Beams (C2) + Cost of pre- 
stressing cables with its sheathing and anchorage in 
transverse Beam (C3) ...3.24 

where for T-section shown in Figure 3*1 

Cl =( a*L + 5-¥^[¥IF(X5-5) - (A - 5-X4)]'\cC-N ...3.25 

and A = area of concrete in longitudinal beam 
¥^ = thickness of transverse web 
CC = cost of concrete/cm^ . 
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02 = [(ON - 1)-CHT¥ + /.DC+CT¥](1 + -|~)L-N + CAC-CN-N 

. . .3.26 

in this ON is total number of sheathing required for cables. 

CHT¥ cost of cable consisting of 12 numbers of 7ra.m dia wires 
with sheathing per cm. 

ADC 008 ^ ' of wires ofiVm.m dia required in the last cable/cm. 

GaC cost of anchorage for cables of 12 numbers of 7m. m dia 
wires per cable. 

CT¥ cost of .one; sheathing per cm 


and C3 = TCN-2b*CHT¥ + TON* CAC ...3.27 

TON = total number of cables required in tranavferse 
direction. 

For a Hollow Box section shown in Fig. 3.2, the different 
costs are 

01 = [A-L + 5.0(X4 - X2 - X1).(WLF - 2.X3)-¥^]CC-N 

. . .3.28 

02 = [(CN-1)-CHT¥ + ADC-CT¥][l + CAC-CN-N 

...3.29 

and 03 = TON- 2b- CHT¥ + TON- GAG ...3.30 


3.5.3 Constraints; The constraints are described here in 
the normalized form, 

3. 5.3.1 Constraints Common to the Two Problems 
(i) Minimum section modulus constraints 


(a) For top fibre 

=^|EiB . 1.0^ 0.0 


...3.31 
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^tmin gi'^en. Toy equation (3.7), is actual section modul 
at top fitre. 


(b) For bottom fibre 

„ ^b min p, , ^ 

^ 2 ”*” 2 < 13 ^" 0 • 0 


...3.32 


^bmin given by equation (3.6), is actual section modulus 
at bottom fibre. 

(ii) Permissible stress constraints for longitudina.1 beams. 

(a) Denoting actual bottom fibre stress at transfer 


condition as f, .. 
f bar 

g = - a.O ^0.0 

^ ^bi 


.3.33 


where f^^^ is allowable stress in bottom fibre at 
transfer condition. 

(b) Denoting actual top fibre stress at transfer 


condition as f , 


~ f4. . 

tai 


1.0 < 0.0 


. ..3.34 


where f . . is allowable stress in top fibre at 


transfer condition. 

(c) Denoting actual bottom fibre stress at working 


condition as f-j^^ 

^b ^ 

g D _ 1^0 <- 0.0 ...3.35 

^ ^ba 

where f-, is allowable stress in bottom fibre at 
ba 

working condition. 
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(d) Denoting actual top fibre stress at woi'king condition 


as f 


taV 


S6 


t a 


- 1.0 


0.0 


. ..3.56 


where is allowable stress in top fibre at working 
condition. 


(iii) Ultimate strength constraint 
M 

grj = - 1.0 ^ 0.0 


M 


...5.57 


us 


where M-,, =2,5M + 1.5M, and M is given by 

11 XX Pu CL XLS 

equation (5.20). 

( iv) Permissible stress constraints for transverse beams 

(a) Denoting the ,a.ctual stress at transfer condition in 
bottom fibre of transverse beam as f. 


'btai 


^8 = 


' btai 

■^bi 


- 1.0 0.0 


... 5.38 


(b) Denoting the actua.1 stress at working condition in 


top fibre of transverse beam as f , , 

gg = _|X&_ _ 1.0 ^ 0.0 


...5.39 


3. 5. 5. 2 Additiona.l Constraints Imposed on T-section 


(i) For under reinforced design 


SlO 


■^"tsw 


- 1.0 ^ 0.0 


... 3 . 40 


0.24 b' d Pc 

jfcermB appearing in equation (3.40) are defined in 
section 3.4. In case of few T-sections 4+^,, given by 
equation (5.18a) becomes negative. As such sections 
are very much under reinforced, this constraint is 
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redefined by assigning it a value equal to 0.0. 
(ii) Constraint on depth of the section 


15iX5 


- 1.0 -s 


0.0 


. . .3.41 


Sll - 1 

where 1 is span of bridge and X5 is overall depth of 
T-section. This constraint restricts the depth to L/15. 


(iii) Heavy top design of T-section 


’12 


- 1.0 ^ 0.0 


. . .3 . 41a 


This constraint will give section modulus in top 
fibre mdre than bottom fibre. 


(iv) Geometric constraints 

13.0 
XI 

14.0 


®13 " 


§14 

§15 


X2 

20.0 

X3 


- 1.0 ^ 0.0 

— 1.0 0.0 

- 1.0 ^ 0.0 


...3.42 
.. .3.43 
...3.44 


These constraints mean that the thickness of top flange is 
to be more than 13.0cm, the thickness of web is to .be- jaore 
than 14.0cm and thickness of bottom flange is to bermorh than 
20.0cm. 

Thus there are a total number of fifteen constraints 
in T-section. 


3. 5. 3. 3 Additional Constraints Imposed on Hollow Box Section 


(i) Bor under-reinforced design 


^10 " 


A, Bs 
tsw 


0.24 b ' d Bs 


- 1.0 :< 0.0 


... 3. 45 


The terms coming above are defined in section 3.4. 
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(ii) Constraint on overall depth 

15iM 


- 1.0 


0.0 


^11 ” I 

This constraint restricts the depth to L/15. 
(iii) Geometric Constraints 


. . 3. 46 


. _ 14.0 

= 12 - XI 

r _ 14.0 

=13 " X2 

14.0 

X5 


Sl4 


- 1.0 ^ 0,0 

- 1.0 -^ 0.0 

- 1.0 ^ 0.0 


..3.47 

..3.48 

..3.49 


The thicknesses of top flange, bottom flange and web 
is to be more than 14 cms. Thus there are a total of 
fourttaen constraints in the Hollow Box section. 


3.5.4 Formulation of Optimum Design Problem 

Thus the optimiom design problems for T-section and 
Hollow Box section bridges turns out as; 

(i) T-Beam Section 

Minimize 0BJ(X1, X2, X3, X4, X5) 
subject to 

g. ^0.0 ; j = 1, 2, ..., 15 

(ii) Hollow Box section 
Minimize 0BJ(X1, X2, X3, X4) 
subject to 

g^ ^ 0.0 j j = 1, 2, . . . , 14 

where objective function, OBJ and constraints g., for each 
case are defined in sections immediately proceeding. 
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3. 5,5 Parametric Study 

Parametric study has heen done to have an idea of change 
in cost. The parameter for hoth the problems is number of 
longitudinal beams. This number changes f rom 5 to 9 for 
both the sections. It is suggested [7] that the number of 
beam has to be more than 4 for the orthotropic plats analysis to be 
justified. That is why the minimum number of beams considered 
is 5. 


3.6 Ontimization Technique 

The interior penalty function approach [5] is employed in 
solving the non-linear programming problem of finding the 
minimum of OBJ given by equation (3.24) subject to constraints 
g. ; 3=1, 2, M, given by equations (3.31) to (3.49). 

J 

M is the total number of constraints. As there is no equality 
constraint to be handled the penalty function is defined as 


9 (X, r) = OBJ - r *2: ^ 

3=1 Sj 


...3.50 


where r is penalty parameter. 

Sequential unconstrained minimization technique is used 
to find the minimum of ^(X, r). This essentially means 
finding minimum of 9 (X, r) as r sequentially goes to zero, 
and in such a situation minimum of 9 (X, r) is minimum of OBJ, 
The imcons trained minimization of 9 (X, r) has been 
pefformed by the algorithm due to Davidon Flecther P/owell' 
[22], In this method the move is generated as 
\ = -[H]^ M. 


...3.51 
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Df^ is gradient at point ; [h] is a positive definite 
matrix of n x n size, where n..,.is number of variable in the 
problem. To start with [h] is an identity matrix. 


The new point X^ _ is obtained by 


"^i+1 


X. + 6 . • S. 
1 *^11 


where is a scalar greater than zero chosen to minimize 9 
function along starting from X^. The matrix [h] is 


updated in each iteration as 


BJi.i 


^ ( B . • S . ) • Y . 

1 1 


[H].(Y. )(Y^)^[ h]^ 

(Y^)'^[H]j_(\) 


...3.55 


where Y. =Df. . -‘Df. ...3.54 

1 1+1 1 

The linear search performed in to obtain was 

done by Fibonacci method [2l]. This makes use of the 

sequence of Fibonacci nuinbers for placing the experiments in 

the initial range of uncertaincy, and reduces the range by 

0.00103 with help of fifteen experiments. The sequence of 

Fibonacci numbers are given as 

F = F-, = 1 
0 1 

and F^ = ^n-1 ^n-2 » n = 2, 3, 4, ... ...3.55 

yielding thd sequence 1 , 1 , 2 , 3 , 5 , 8 , 13, 21 , 34, .... 



CHAPTER 17 

RESULTS AND DISCUSSIONS 

4 . 1 Introduction 

The minimum cost design of two bridge decks each having 
different types of cross-sections of the longitudinal beams 
has been carried out in the present work. The illustrative 
examples of spans 24 ffi and 30m in each case have been presented. 
As a result of the present work a. computer program has been 
developed for minimum cost design of bridge decks. The 
listing of the prograni and the subroutines in fortran 17 is 
given in the appendix. 

The computer sub-program referred to as TEEBRI performs 
the minimum cost design of bridge decks constisting of a 
number of T-section longitudinal beams, while the computer 
program referred to as BOXG-IR is for minimum cost design of 
bridge decks consisting of a number of Hollow Box sections. 

The formulation of these problems is given in Chapter III. 

Cost of different items required to find the cost of 
bridge deck is given in Table 4.1. The cost table has been 
taken from reference 9. The allowable stresses in concrete 
and steel wire is presented in Table 4.2. This table 
has been generated from reference 5.^ 
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Table 4.1 : Cost of Items 


Item 

' . , i 

Descriptiorx | 

Rate 

■" ' "■ ■ ’ 

Concrete 

1 

I 

i 

\ 


(a) 

j 

Grade of concrete M-350 | 

265.00 Rs./m3 

(b) 

Grade of concrete M-425 

297.00 Rs./m^ 

H.T.S. 

i 


wires 



(a) 

Gable consisting of 12 nos. of -7m.m | 

( 



dia wires of H.T.S. with sheathing 

[ 

1 18.00 Rs./m 

(1) 

One un tensioned wire of 7m. m dia 

0.91 Rs./m 

(c) 

Two anchorage for cable consisting 



of 12 nos. of Tm.m'dia wires 

110.00 Rs./cabl 

(d) 

Cost of sheathing for cables con- 


' 

1 

sisting of 12 nos. of 7m. m dia wiresj 

7.00 Rs./m 
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Table 4.2 ; Allowable Stresses 


Item ! 

i 

- f 

Description of stress | 

i 

Nota- 

tion 

Permissible | 
Values 1 

Concrete | 

1 

i 

1 

i 

3 


1 

! 

I 1 

(a) j 

Allowable stress in bottom | 

^bi 

j \ 

i 

1 

fibre at transfer j 



i’ 

(i) For concrete grade M-350 | 


1 140.0 Kg/cm^ 

1 

i 

(ii) For concrete grade M-425 | 


157.14 Kg/cm«j 

(b) 

} 

Allowable stress in top fibre | 
at transfer ! 
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(i) For concrete grade M-350 j 


0.0 Kg/ cm^ 


(ii) For concrete grade M-425 1 


0.0 Kg/ cm 2 
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1 Allowable stress in top fibre | 
at working j 

(i) For concrete grade M-350 j 
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110 .0 Kg/ cm^ 


(ii) For concrete grade M-425 i 
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Allowable stress in bottom i 

fibre at working 1 

‘ (i) For concrete grade M-350 
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0.0 Kg/cm2 j 
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, (ii) For concrete grade M-425 


0^0 Kg/ cm 2 1 
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Steel 

1 Ultimate strength of steel 
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}: wire 
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1 i 

! i 

15000 Kg/cm2 1 
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4 . 2 Corn-par is on of Results 

The sections and spans taken up for study in present work 
also appear in the list of standard sections recommended by 
SSRG handbook [9]. Furthermore, the sections corresponding to 
minimum cost of bridge decks have been recommended in the 
said handbook. The optimization technique used in finding 
these optimum sections in rather crude. In the referred work, 
the cost of bridge decks with beams increasing in depth by 
one centimeter is compared and the depth which corresponds to 
minimimi cost is recommended as optimum design. This essentially 
is formulating the optimum design problem r.with one variable. 
Moreover the cost of transverse prestressing cables has been 
neglected. As compared to this the computer programs 
developed in the present work handles the same optimization 
problem having five variables and fifteen constraints as 
fomulated in Chapter III. In order to compare the results 
of the present work with that documented in reference 9 the 
sub-program TEBERI was modified to solve the same problem but 
having five variables and thirteen constraints. The two 
constraints dropped are gg and g^ which go to check the 
stresses in transverse beams as given by equations (3.38) and 
(3»39). The objective function was also modified to exclude 
the cost of transverse prestressing cables. Ths results 
given in Table 4.3 are for' 24m span bridge deck of concrete 
grade M-425, the number of longitudinal beams considered are 



Table 4,’3 ; Comperlson of optimxim design v/s standard section design of SSiiC 
(Span 24m, Concrete grade M-425) 
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six and seven. It can be concluded from Table 4.3 that the 
sections given in reference 9 do not q^uite correspond to 
minimum cost design. The present work has shown the reduction 
in cost by 8.4 percent in one case and 9.15 percent in other 
as is seen from Table 4.3. 

4. 3 Results of Parametric Study 

The results of the parametric study obtained for the 
two illustrative examples as formulated in Chapter III can be 
categorized into the following. 

(i) 30m span bridge decks comprising of T-section longi- 
tudinal beams, number of which changes from 5 to 9. 
(Table 4.4, seria.! number 1 to 5). 

(ii) 24.m span bridge decks comprising of T-section longi- 
tudinal beams number of which changes from 5 to 9. 

(Table 4.4 serial number 6 to 10 ). 

(iii) 30m span bridge decks comprising of Hollow Box section 

longitudinal beams, number of which changes from 5 to 9. 
(Table 4.5, serial nmber 1 to 5). 

(iv) 24m span bridge decks comprising of Hollow Box section 
longitudinal beams, number of which changes from 5 to 9 
(Table 4.5, serial number 6 to 10 ) . 

In the above results the grade of concrete is M-350. 
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Table 4. 5 : Hollow Box Section Bridges 
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4* 4 Discussions of Results 
(a) T“Section Bridge; 

In general it can be said that the cost of bridge increases 
with increase in number of longitudinal beams. Let us consider 
the 50m span bridge. The cost of 5 beam bridge is minimum. 

This cost is 0.926 times the cost of 7 beam bridge. The 
weight of one beam of 5 beam bridge deck is about 1.22 times 
the weight of one beam of 7 beam bridge deck. When the 5 beam 
bridge deck is compared with 9 beam bridge deck its cost is 
0.8 times the cost of 9 beam bridge, and the weight is about 
1.29 times the weight. 

Taking the results of the 24ffi span bridge the cost of 
5 beam bridge deck is about same as 7 beam bridge deck. But 
if weight of individual beams are compared the 5 beam bridge 
is about 1.58 times the weight of 7 beam bridge. Again 

comparing the 5 beam, bridge with 9 beam bridge it is 

observed that 5 beam bridge is about 0.76 times the cost of 

9 beam bridge and in weight they are 1.18 times heavier. 

The weight of individua,! beams can change overall economy at 
the time of launching of the bridge. A heavy beam will need 
more careful handling. Thus if in the span under consideration 
(24m and 50m) a 7 beam bridge deck is constructed there will 
be an average increase in the cost of bridge by about 4 percent, 
but the weight of individual beams to be handled will be 
about 2'5 percent less. Since the detailed data regarding 
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handling and launching of the heams was not aTailahle a 
rigorous study could not be pursued, however, it is felt 
that the 7 beam bridge deck will turn out to be over all 
economical. 

The active nature of the constraint g 2 means that the 
optimum design point corresponds to a heavj?' top flange design 
with the section modulus at bottom fibre held near the 
minimum section modulus required. This will keep the minimum 
eccentricity of the longitudinal pre stressing force given by 
equation (3.9) within the depth of the cross-section. 

The active nature of constraints eaid naean that at 
optimum point, two out of the four stress inequalities given 
by equation (3.1) to (3.4) are critically satisfied. Cons- 
traint checks the bottom fibre stress at transfer, whilfe 
checks the bottom fibre stress at working conditions. 
Constraint g^ is active for all designs as the longitudinal 
prestre'ssing force has been calculated by assuming that the 
bottom fibre stress at working condition has reached the 
allowable stress. 

The constraints due to minimum thickness of web and the 
maximum depth of the section are also active. As a combined 
effect of the two, the tendency of the variables at optimum 
is to decrease the web thickness and increase the depth. 

This is because the moment of inertia of such sections are 


higher. 
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(b) Hollow Box Section; 

In case of Hollow Box section bridges the cost of bridge 
deck increases with increase in ninaber of longitudinal beams. 
In the case of 50m span, 5 beam Hollow Box bridge is about 
0.792 times the cost of 9 beam bridge deck. Bor 24m span 
brid^^e this ratio is about 0.845. Thus in general for the 
Spans under study 5 beam Hollow Box bridge is economical. 

When the cost of Hollow Box bridge is compared with the 
equivalent T-section bridge in number of longitudinal beans, 
then the Hollow Box bridge turns out to be uneconomical. In 
30m span 5 bCcam T-section bridge is about 0.875 times the 
cost of 5 beam Hollow Box section bridge. This ratio is 
0.818 for 24m span. 

However in certain situations like minimum head room 
consideration, or when an old bridge is to be replaced and 
the expence regarding existing road leading to the bridge 
is to be avoided, Hollow Box section has to be preferred. 

This is because the overall depth of cross section is less 
in Hollow Bax section when compared with the overall depth of 
T-section bridges, for example, in case of 30m span bridge 
the depth of 5 beam Hollow Box section bridge is 0.91 times 
the depth of 5 beam T-section bridge. This ratio for 2 4m 
span bridge is 0.725. 

Again the constraint g^ which checks the bottom fibre 
stress at working is always active. This once again is 
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because the pre-stressing force in longitudinal beans has 
been calculated assuming a section with the bottom stress at 
working to be critically satisfied. The optimality of the 
Hollow Box section is governed by the minimum thickness 
constraints. The minimum thicknesses of top flange, bottom 
flange and web have been prescribed to be above 14 cms. 
constraints due to them are active at optimum. 


4.5 F ew Observations on Numerical Aspects of Optimization 
Routine 

The optimization of penalty function is performed for a 
successively reducing value of penalty parameter r in SDMT. 

A large value of r will require too m.any unconstrained 
minimization cycles. As suggested by Fox [20 J the starting 
value of r can be taken as 

r . 


M _1 
j=l ^ 


where OBJ is objective function, g. is constraint and 

J 

M is total number of constraints. They are as defined in 
Chapter III. In the present work the value of r given by 
equation 4-1 was very high. By several trials it was found 
that a starting value of r = 80 with a reduction factor of 
0.2 to find the subsequent r, gave satisfactory results. With 
these values the unconstrained minimization cycles required 
was four to five. The results given in Table 4.6 are from 
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computer output when TEBBRI was optimizing 50m span Bridge 
deck donsisting of 5 numbers of longitudinal beams. In this 
only four unconstrained minimization cycles were required to 
arrive at the optimum point. 

Alongwith a reasonable starting value of r, a proper 
convergence criterion is also required to identify the 
optimum point and terminate the iterative scheme. T 

The D.F.P. algorithm used for unconstrained minimization 
is terminated whenever any one of the following conditions 
is satisfied. 

(i) Df • S ; where Df is gradient vector and S is the 

move vector, 

/ i 1 

(ii) — = ; where k is iteration 

r) 

! - 1 . number of linear search. 


( iii) max | DX| ^ ; where DX = ^ j 

X^ and are optimimn point arrived at after k"^ 

k-l"^^ linear search. 

(iv) I)f^-Df<e. 


The SUMT algorithm is terminated when both of the 


following conditions are satisfied: 

<p(X r ) - 0BJ(X. ) 

( i ) I__i i 1_ 

OBJ(X^) 

where X^ and are the value of optimum point and penalty 

ii 

parameter in the i -unconstrained minimization. 
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2 ^ are 

th th 

optiraxmi point arrived at after i and i-1 uncons- 
trained minimization cycles respectively. 

The value of the small positive numbers taken in the 
present work are as follows: 

= 0.0005, eg = 0.0001, = 0.0005, e^ = 0.0001, 

= 0.001 and = 0.001. The effect of these convergence 
criterion are presented in tables 4.6 and 4.7. The i-'results 
given in Table 4.7aro for the same problem for which the 
results are shown in Table 4.6. 

From Table 4.7 it can be seen that the reduction of 
penalty function is maximum in the fourth linear search. In 
the fifth linear search the reduction in penalty function 
is very low and hence D.F.P. algorithm is terminated. The 
results of Table 4.6 shows the change of variables to 
optimum with successively reducing values of r. After end of 
fourth -unconstrained minimization there is no change in 
design vector frbh the design vector of 33^d iteration and 
difference between penalty function and objective function 
is also very small, hence the program declares the current 
point as the optimum point. 


(ii) max |DXl~-'e^ : where DX = X. - X. , : X. and X. 

O 1 1—1 1 X 




TaVle 4.6 ; Flow of Design Variables for Different 



















Table 4.7 ; Change of design variable in different linear search (r 
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4. 6 C onclusions and Recommendations 

The cost of bridge deck consisting of either T-sections 
or Hollow Box sections increases as the nviraher of longitudinal 
beams increase. T-section bridge is always cheaper r of the 
two sections and hence its use is recommended. Even though 
the cost of 7 beam bridge is about 4 percent higher in 
comparison with the 5 beam bridge deck, use of 7 beam T-section 
bridge deck is recommended as the modules are lighter by 
about 23 percent. However, the use of Hollow Box section is 
recommended when the overall depth of cross section has to be 
on lower side. 

At the time of formulating the optimization problem for 
minimum cost design of bridge deck the variables have to be 
properly identified along with right choice of objective 
function. By over simplification of the formulation of the 
bridge deck problem and neglecting the cost of cables required 
for transverse prestressing the sections recommended by the 
Handbook of S.B.R.C. [9] do not correspond to the minimum 
cost design. 

4. 7 Scope for Future Study 

In the present work the local effect of wheel load in 
the top slab has not been accounted for. This, being a 
secondary stress distributioh, it can be superposed on the 
stress distribution produced in the bridge as a whole to 
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give the resultant stress distribution. There are several 
methods available, in the Ixtefettorb' \ to study these 
effects. They can be included to check the optimality of 
this problem,. 

The minimum cost design of these sections have to be 
taken up for a wider range of sp'an to have a clear picture 
of obehaviour of optimTim point. Other commonly used sections 
for bridge deck i.e, inverted T-beam (low span range) and 
composite I-section could also be studied for standardising 
the minimum cost design of bridge decks. 
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c. 

C THIS TMCrRPORATES THE SEQUENTIAL UMCOMSTRAINFD t xjT :^r Z A •"■ HU 

c technique 

C SPAN=SPAM OF THE BRIDGE 

C MB= number of LONGITUDINAL BEAMS IM THE BPTDGF DECK 

C WLF = WIDTH OF TOP FLANGE OF LONGITUDINAL FT AMS 

C WTF* WIDTH OF TOP FLANGE OF TRANSVERSE BEAMS 

C THETA IS THE FLEXURAL PARAMETER 

C TPL IS ARRAY STOPTNG TRANSVERSE POSITION OF LOAD FO^ 

C TRANSVERSF RFNDJ NG MOMENT 

C PL IS ARRAY STORING POSITION OF LOAD IN LONG! tuOI NAI. D^RFST'O! 

C FOR MAXM. TRANSVERSE 8.M. 

C WP70R IS ARRAY STORING WEIGHTING FACTORS FOP Tjp LOALOMC 

C WPFP IS ARRAY STORING WEIGHTING FACTORS FOR CLASS AA LOA''^^ MG 

C B IS HALF WIDTH DF BRIDGE DECK 

C 8MFLL IS MAXM.B.M. DUE TO LIVE LOAD ON FOOT-PATH AS PEP ' RC COOr 

C 8MFDL IS MAXM, DUE TO DEAD LOAD OF FOOT-PATH 

C BMWC IS B,M. DUE TO WEARING COARSE 

C BMLL IS B«M, DUE TO 70R/ CLASS-AA LOAD IN BRIDGE SPAN 

C FPTTP IS ALLOWABLE TOP STRESS IN CONCRETE AT TRANSFER 

C FPTWP IS ALLOWABLE TOP STRESS IM CONCRETE AT WORKING 

C FPBTP IS ALLOWABLE BOTTOM STRFSS IN CONCRETE AT WORKING 

C FPBTPIS ALLOWABLE BOTTOM STRESS TN CONCRETE AT TRANSFER - 

C SCU IS ALLOWABLE ' , 

C SCU IS 28 DAYS CRUSHING STRENGTH OF CONCRETE 

•C CHTW IS COST OF CABLE WITH 12 NOSe OF 7MM WIRES WITH SHEATHING 

C PER CM. RUN 

C CC IS COST OF CONCRETE PER CUBIC CENTIMETRE ■ , 

C CAC IS COST CF ONE SET OF ANCHORAGE PER CARLE , , ; 

C ■ CITW IS COST OF ONE WIRE OF 7MM DIA, 

C CS IS COST OF SHEATHING PER CM 
C RE IS LOSS IN PRESTRESSING FORCE 

C X IS ARRAY STORING DESIGN VARIABLE 

G • S IS ARRAY STORING THE CONSTRAINTS 

,C M IS NUMBER OF CONSTRAINTS , N IS NUMBER OF VARIABLES ; 

C 

C , - ' ■ • 

:OMMnN/SETl/SPAN,NTF, WLFf NBf THETAtT’Lf PL, P, WPTOR » HPFP, 8, BMFLL, BMFO 
1 L,BMWC,BMLL, FPTWP, FPTTP, FPBTP,FPBWP,, SCU, FS,CC, CHJW, CAt,G'!<S,RE 
?,C1TW,CS 

. C0MM0N/SET3/PT,P0,NTC,NCfAPC. . ■ ^ 

DIMENSION X( 10) ,XDf lO),Gll51iTPLt'4I, PLC'7),P<TI,WPT0RtP'),WPFPC D) 

. I E P SM { 4 ) , XS U 0 ) , ' ' ■ ^ 'T"'"- " . , : ^ 

DATA TPL/l5B.0,79.0,0*0*-79»a/ ' ‘ /-.i-V"’ 

DATA PL/302«.0, 698.0, 850.0, 1063. Q,1230,0,i5tf5.0, 1,642.0/' ' ' 

'^'ATA P/8000.0, 12000.0,12000.0, IT'OOO.'Orl'TOOO. 0, ifOOO'.O^ ITOOO. 0/' , . 

■ ■ DATA WP70R/1.4l85,1.783.Z,'0.?194,'0.0,0»0,,<J#07t9,0«TH,0*0',0,'0/. r'.'' , 



o o o o o o 




HATA WPFP/0.0,0.0,0«0417,0.'^&33,0.375,O.D,'-u:iil7f 
OATAIMB,ND,B,PE,FSfCHTW,CAC,ClTWt CS/f ,5,55 1"^- < ■. >. I : - } ] 

i.0,n«')09i,noU7/ 

data FPTWP,FPTTP,FPBTp, FPBWP,SCU»CX/110oO,OaO ’ s : i-.,-* U 
10.000265/ 

DATA SPAM,R^*LL^WC» FDL/ 2400 .0, 0,43833 5E 08, 134b-), ^91.3/ 

DO 5 T=lt4 

5 TPLny=TPLn )/b 
WTF=SPAN/4.0 

FLL=400.0-( C 40, fj*SPAN) / 100,0-300, n)/9,0 

FLL=FLL/10CtOruO 

BMFLL=FlL*SPAN=^'SPA^!/4.0 

BMFDL=FDL417 5, Oi««SPAUt.SP AM/40000.0 

PKWC=WC=?'750«0*SPAM*J PAN/80010.0 

data CiIMAXi Ef N/0»2,lO,D,5fc-"02f 15, 5/ 

M=14,N=4 for hollow BOX GIRDER 

DATA X/20«Ojl6»D,40»0,80.0,149,G,D,3fDBU,OBOfUtO,(*i»0/ 

FOR BOX GJRDEP X U ) = 1 S. 0, X ( 2 ) = 16. 0, X ( 3 ) = 1 5 . 0 , X { 4 ) = 1 " 0« 0 i 


FOR COST DESIGM OF PRF- 

STARTIMG P3IMT IS*) 


STRESSFD 


DO 3 1=1, N 

, 3 xs{n=x(i) 

7 PRINT 500 

500 F0RHAT{5X,*stART OF SUMT 
ICONC. BRIDGE*,/, 20X, *THE 
WLF=2.0*B/FLOATCNB) 

DO 8 1=1, N 

8 xn)=xs(j) 

PRINT 510,{I,X(I ),I=1,N) 

ITER=0 

PRINT 450,SPAN,NP,SCU,FPTWP,FPTTPf FPBTPfFPBWP 
450 FORMATt 10X,*SPAN IS*,F10. I, 2X, *CMS*, /, lOX, *NQS,OF BEAH IS*, 13,/ 
110X,*CUBE STRENGTH OF COMC. IS*,Fa,l,* KG/CM-S3*, /, $X, ^ALLOWABLE - S' 
2TRESS TOP{WORK|NG)*,F10.1,/ ,5X, *ALL OW ABLE STRESS TOP ( TRANSFER ) *, 
4F10. 1,/,5X,*ALL0WA8LE STRESS BOT. (TRANSFER 1*, F 10. 1, /, 5X, *ALL0WAB'LE 
5 STRESS BOT. CWQRK!NGJ*,Fl0.iI ' '■ ' 

11 R = 0. 0 

CALL FTNCG,X,FI ,nEJ,«,N,R> 

R = A8SC0BJ/GKS) 

IF(R.GT.BO) B = 80.0 
PRINT 512,R 

510 FnRMAT(5(5X,*X(*,I2,*J =*,E16,8/M 
512 FORMAT(////,50X, *8=*, 616.617 


S', 


ili 


PRINT 513,IC{I>,J=l,M> 

PRINT 515,PT,Pn,NTC,NC, ARC 
515 FORMATIZE20.B,2IIO,F2D.8) 

IFCFI,GT,0.5E 09) GO TO lOil 
00 10 I=1,N 







o o o o o 


13 xna )=x(T ) 

Fl=OBJ-P=^GKS -♦ 

PRINT nSJ 

514 F0RMAT(5X,*F =*» £16*, P, 5X, *OBJ =*,F1S.8) 

25 TTER=ITER+i 

CALL UNCHNf N,M,X»G» Ef TMAXt R) 
call FTN{ GtXfFf OBJf M, N,R) 

I F{F,GT.P®‘^F 09) Gn TO 1001 
F1=:ABS(F-F1)/aBS(F) 

F2=0«0 

DO 30 J=1,M 

A=XD{I)-X(I) 

IFCABSf A)»GT„F2} F2=A 
30 CONTINUE 

r F{E1.LT.O.O(.3 . ANn,~2.LT.0.nOl) 30 TO :;0 

PRINT 515,PTsPQsNTC?NCj APC 

RsC*R 

PRINT 512, R 

JFflTER.GT.lO) GO TO 60 
DO 40 1=1, N 
40 XDd )=X{n 
F1=F 


50 

605 

65 


60 

606 

513 

1001 


1011 


100 


GO TO 25 
PRINT 605 

F0RMAT(5X,*CnNVERGENCE ACHIEVED*) 

PRINT 510, (I, X(! ),I=1,N) 

PRINT 513, (G{J), I=1,K) 

PRINT 514, F, OBJ 
SO TO 1001 
PRINT 606 

format(5x,*convergenc£ not achieved*) 

GO TO 65 

F0RMAT(5X,*THE constraints ARE*,/9E12.4,/6E 12-4) 
MB=NB+1 

IFCNB.GT-9) GO TO 100 

GO TO 7 

X(3)=X(3)“2*0 

IF(X{3).LT,20.0) XI3)=20.0 

X(4)=X(4)*2.0 
TK=I K+1 

TF(TK.GT.30) GO TO 60 : 

GO TO 11 

STOP •- . ' 

END ' C;:' 

SUBROUTINE FTN{G,X, F,QBJ, H,N,R) . 


::4. 


, , ' 


iiHli: 








* -'i !• ' ' ' 


THIS EVALUATES THE PENALTY, FUMCTI0N( iT,OBjeCTlVe.''FUNGTl ONC OBJ) 
CONSTRAINTS ^ AT A .POINT ,X TOR; A'^tTEB 7: ' 

* "■ ' ' ■■ — — ^ 



o o o o o n 


MOTE THIS SUBROIJT^ME IS FOR T SECTIOM 
SUBROUTINES PEOUlRFO ARE MAIN, TORSIO ,GEN, 01 STr 
AND FJBfiO 


■p ASP" , U>\,CD'i 


COHMON/SeTl/SPANtWTF,WLF»NB,THETA,TPLf PL f P f P'PTDR ^ P P, 
lLtBMWC,B“^LLtFPTWP,FPTTP,FPeTP,FPBWP, SCU» FS , CC? PH^Hf CA! 
2,C1TW,CS 

:0MH0^4/SET3/PT, PC, NTC,NC» APC 


P? ‘-’^'FLL. '"-"F' 
f [7 K, , R E 


01 MENS! QN X( 10) jGf 1 L) , K0{ 5 , 0 ) , KUt 5,9 ),MEU< 3, 4 ) , , 4) , l d7,'R » ^ . 

1WPFP(9) ,0K(9,9), APK0(9),ABKU(9),TPL< 4),EIH6),PH“:"),P( 71 
REAL KO,KU,FFW,MEU 
HU=X(5)’-X(1)-X(3 ) 

APC=WLFTX(1) -J-HUTXf 2)+X(3)TX<4) 

FNPCB=XC4)>!'X(3)>^X{3)TD«5+X(21=?'HJ=5(X( ?) +01)40, f> ) + x { ;■ ):''WLF*(X( M-Xf 1 ) 
1*0.5) 

YPB=Ff'*PCB/APC 

YPT=XC5)-YPR 

CHI = WLF*f XU )**3)/l ?.0+WLF*X(l )*( { YRT-X{ 1 ) + CHU4*: UUXf 2) / I 

12. 0+HU* X { ? ) * { ( YPT-X C 1 ) -HUTO. 5 ) **2 1 +K C A- )F=C X ( 3 ) ** 3 ) / 1 7 , 0+X f 4 1 * Xf ’ ) T < 
?{ YPB-XUl^O.SlXc^Z) 

ZPB=CMI/YPB 
ZPT=CMI/YPT 
S(121=ZPB/ZPT-1.0 
RATIC=WLF/X( 1) 

CALL TORSln( RATJ 0) ' r;;:; ■ ■ 

OI=RATIQ*WLF*{X( 11**31/2.0 
RATIO=HU/X{ 2 ) 

CALL TORSIOCRATIOl . . 

OI=OI+RATIO*HU*{X( 21**3) 

RATI0=X{4)/X{3) 

CALL T0R5I0 CRATIOl U' 

0I=01+RATI0*XK)*(XC 31**3) 

WTB=25.0 
RATIO=WTF/X(l) 

CALL TORSIOI RATIO) 

Oj=rATIO*WTF*{X( l>**3)/2.0 t 

HI=XI5)-XU)~5.0 
RATIO=HI/WTB 
CALL TORSIO(PATIO) 

DJ=04+RATin*HI*t WTB**3) ,, 

ATB=WTF*XU ) +WTB*H| ■ 

T M8= WT8 * { HI **2 1 *0. 5+WTF*X {„l J *{.Hl +X U 1 *0, 

3TB=thr/atb -'U . 




k. 


0TT=Xt5}-DTB -5.0 






TJ=WTB*(HI**3)/lZi0+WTB*HI*l,CbT6-Mr*0,f,1*^;'2>+M7F*lXC 1)**3)/12. 4+WT 


1F*XU)*( (DTT-XU )*0.51**2> 
ZT8B=TJ/0T8 
ZTBT*TJ/DTT 






A ^ "it ' , 

-P X'M-- 

“ 7 v7> ’ ‘ 





65 


SI=CMT/WLF 

S|0=0ir/WLF 

J=TJ/WTF 

SJO=OJ/WTF 

THETA=R*C CSI/SJ)=<'«0«?5)/SPAM 
ALPHA=D«2175’i‘(Si:0+Sj0) /SORT{ Sl^i'SJ ) 
SALPH=SOPT( ALPHA ) 

CALL niSTRT (THFTA, ALPHA,KO,KU} 

CALL GeN(KOtCK) 

DO 10 I=sl»9 

XX=0,0 

DO 5 J=l,9 

5 XX=XX+WP70R( J3*0Kf l5 J5 
10 ABKO(I) =XX/4eO 
CALL GEN(KU,0K3 
DO 20 T=i, 9 
xx=D.a 
DO 15 J=l,9 

15 XX=XX+WP70R( J)*OK(I, J) 

20 ABKUd) =XX/4.0 
Y=0. 0 

DO 25 1=1,9 
XX=ABKU(I)-ABKn{I) 

ABKO(I)=ABKO{J ) + XX’«'SALPH 
IF(Y.GT*ABKO{T ) ) GO TO 25 
Y = ABKO{I) 

IJ=I 

25 CONTINUE 
I M=I J-1 

0IFB=FL0AT{IM)*E/4«D 

XX=WLF/2.0 

MS=NB-1 

DO 30 1=1, NS 

NBS=I 

Y=OI FB-XX 

IF(Y,LT.WLF) GO TO 40 
30 XX=XX+WLF 

40 DKVJ=ABK0{IH)+4.0*Y*( ABKOC IJI-ABKO? IM >)/B 
CALL GENCKOjOK) 

DO 45 1 = 1,9 
XX=0,0 
DO 42 J=l,9 

42 XX=XX+WPFP( J)40K(I, J) 

45 ABK0{iy=XX/2»0 
CALL GEN {KU,OK) 

DO 50 '1 = 1,9 _ 

X X=0 • 0 ■ ' - ■'■'>'■1 ‘ 

DO 48 J=l,9 

48 XX=XX+WPFP( J340KM, J> 

50 ABKU(I)=XX/2.0- ... 



DO 55 1=1,9 
XX=ABKU(I )-ARKO( n 
55 APKOC I) =ABKO(7 )+XX^SALPH 

DKF=ABKn( JK)+A«0«( ABKndJj-ABKOnx )) /B 
DWD=(WLF*{X t 5)-5,0)“{ APC--X{4)*5«0} 002^ 

WS = APC«'0.00 24 

BN1S=WS*SPAN*.SPAN/8.0+-DWD*SPAN/2.0 

BM=FLOAT(NB) 

f e LL= 1 . 1 *0 KF MFL L/ 8N 

PB»^FDL=1.1*0KF=XBMFDL/BN 

PBVWC=BMWC/BN 

PBB'LL=BMLL=«'lolX'l,15^nKW/BN 

peMA=PBMFLL+P8MFDL+pBMWC+PBHLL ■ 

PBMD=BMS+PBHFDL+PBMWC 

pbml=pbmfll+pbmll 

ZTP=(PBMA+{ l,.’}-leP/kF)X=BHS)/( FPTWP-FPTTP/iU ) 
ZBP= SP8MA+{ 1«0-1,0/RF)*BMS) /{FPETP/^E-FPPWP ) 
3(1)=ZTP/ZPT-1,0 
G{?J=ZBP/ZPB”1.'J 
= P-YPB-0.0B5*X(5) 

P0=( PBHA+BMS+FPBWP*ZPB) ^RE/( EP+ZPB/APC ) 
FPBTA=P0/APC+Pn*5P/ZPfi-BMS/ZPB , 

G(3)=FP8TA/FPBTP-1.H 


FPTTA=PO/APC-Pn«ep/ZPT+BMS/ZpT 

5(4)=FPTTP/FPTTA-1.0 

FPBWA={ PO/APC+Pn*FP/ZPB)/RE -RMS/ ZPB-P R^'A/ Z P8 
G{5)=FPBWP/FPBWA-1.0 

FPTWA={ P0/APC~P0*FP/2PB)/RE +BMS/ZPr +PBMA/ZPT 
G (6) =FPTWA/FPTWP-1,0 
CN=P0/<0.6*FS*0.3P5*12.0) 

XX=AMOD(CN,1.0} 

IFCXX.GT.0.0) GO TO 60 
MC=IFIXICN) 

ADC=0,0 
Y=0, 0 
GO TO 65 

y=i, 0-xx 

ADC=XX’}'12,0 

C N»C N+Y ■ ■■ 

MrsTFTXIPMI . 

ATS= {12.0*{CM-Y) +ADC)Xc0.385 - . 'Iv' • 

FFD=YPT+FP ‘ 

UHR=1*5*PBMA+2«5=S'P8KL 


ATS= 112.0*{CM-Y J +ADC)=®=0.386 - • .i- -./{v, -r/* ' ■ - . 

FFD=YPT+FP 

UHR=U5*PBMA+2«5=S'P8KL ^ 

ATSF=0«68*SCU*(WLF-XC2H*XUI/FS ^ •' ' ’ 

ATSW=ATS-ATSF ' ' '' 


G C10)=ATSW#FS/ (0,24«Q*68*SCU*X<2}*EFOJ*-l.a"i:-4i|K4?'r5-r'%=' ;.>■ 

^ ;ti1- , r -i' ■' 


IF ( ATSW*LT.O.O> G { 101 '■ .'P 

USM= ATS W*F S ^^'EF D* ( 1 . 0-0, •rS^'AT SW*»'F S/ 1 X 12 1i^*EFf1<^'SG J 
IFlXd) .GE.EFO/4.0) .GOTO 70 
U S M= US +0 . 7 =«' SC U* X ( I ) « W t F-X I EFO -0*. 5*.X' U H 


•T i ■' V/'.' 
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70 S(7)=UHR/U'^M-1.0 

CALL TRAC0E(TPL,'VEU, mew, THETA) 

DO 80 T=1 ,6 

xx=o»o 

IF{I.GE.4) GO TO 12 
no 59 J=l,4 

59 xx=xx+Mewn ? j) 

EU(I)=XX 
GO TO 80 
72 XX=0.0 
T0=I-3 

no 11 K=l,4 

77 XX=XX+MEU{Tn,K) 

EU(I )=XX 
80 CONTINUE 

DO 82 1=1, 5, 2 

EUd )=EU(T) + f EU{ I + 1)-EU(I ) )*SALPH 
82 EU(I+1)=0, ) 

PI=4.0*ATAN( 1.0) 

YML=0.0 

IK=0 


Si 


DO 90 I =1,5, 2 
IK=IK+l 
E!=FLOATU ) 

XX=0.0 

DO 85 J=l,7 

85 XX=XX+P( J)=!‘S!N(ET*PI*PL( J )/SPAN) 

IJ=IK+1 

YML=EU( I )=^=XX=i‘( C-1.0)=S'=^I J)+YML 
90 CONTINUE 

YMi=YML*2«0*8/SPAN 
TYML=YHL*WTF 

PT=(TYHL+FPBWP*ZTBB+FPTTP*ZTBT/Re)*ATB»=RE/{ ZTBB + ZTBT) 

G{8)=PT/ATB-1»0 

5( 9)= (ZT8T*PT + REX‘TYML)/(RE=I'ATB*ZTBT#FPTWP)-I.0 
PT=5.0=«'PT 

TNC=PT/C0.6X‘FS=«'0«385«12.0) 

XX=AH0DdNC,1.0) 

TFIXX.EO.O.O) go to 100 

Y=l. O-XX - '■ , p 'V 

tnc=tnc+y 

100 NTC=IFIXITNC) 

G(11)=15.0*X{5)/SPAN-1.0 a 
G (13)=18.0/X(l)-l«0 
GC14)=14.a/X{2)-l,0 
s(15)=20.o/xc3)-i.o 

X X* t AP C * S P A N+5 . 0 4w T B’t'f WL F'*t X ( 5 > ^5 . 6 ) -I APC "5 . 0 «« X ('4,1) ))> C C 
XX=XX+CFL0AT(NC)4CS + C ATS*ClTW/12*Or)*SP*H*{U0+8.O*5P*EP/(3. 

14SPAN)) ' ' , ■ ' 

DBJ=XX4BN +TNC*CHTW*2kO*B 


.3*1 


, 'A .. V 

'f ^ 





G^sp an;- 




/ >-C: 

®®sit 
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OB j=obj+tnc*cac+cn*cac*bn 

SKS=0.0 

00 110 1=1, M 

F (G (I ) . EO, 0 ,0) GO TO 110 
GKS=GKS + 1« Vr,(7 ) 

110 CONTINUE 
F = 0BJ-R=^'GKS 
DO 120 1 = 1, 

120 F=0.5E 11 

140 RETURN 
END 

SUBROUTINE tors I Of RAT) 


A: ft 


‘•f’? .'If Jjj-r 


THIS GIVES TORSIONAL CONSTANT FOR A RECTANGUL 


AR SECTTjV 


* * ILi ****-********»*»**....**..***. 

PI=4.0*ATAN{ 1,0) 

X=0. 0 

I F (B,GT» 10, 0 ) GO TO 20 

DO 10 M=1,U,2 
FM=FLOAT(M) 

10 X»X-fTANHCEM^P|*B/ 2,0 ) /(EM=*t5) 

RAT=l,0/3« 0-64, 0*X/n )=«:e) 

GO TO 30 ■ : 

20 RAT=1. 0/3.0 
30 RETURN 
END 

SUBROUTINE TRAC0E( TPL, MEU, MEW, A ) 

THIS generates TWO 3X4 MATRICIES ONE EACH FOR ^0 AWDyUl ■ ^ 

AEU(0*.5} '•'" ^ 

REAL MEW, MEU '' ' - 'iv- 

= = 0.0 

THETA=-A ... 

PI=4,0*ATAN( 1.0) 

B = 1,0 ^ 

DO 10‘T»1,3 
THETA=THETA+2,0*A 

SIG=THETA*PI . : 'v, ...I 

CHG=COSH{SIG) . 



SHG=SJMH{ ST G) 

AM0=PT»THETA/( B’:=S0RT{2,0) 5 
TME=2»a*AMR*R 
SHE=STNH{tme ) 

SNE=STN{T^iP) 

DO 10 J=l,3 
Y=-TPLCJ) 

CHBPY=COSH( AYB^C B+Y) ) 

CS8PY=C0S{ AMB’-M R+y) ) 

SNBPY=SIN( AMB^( B+Y) 5 
SHRPY=STNH(AYB*C R+Y) ) 

CS8PE=C0SCAYB*{B+E) ) 

CHBME=CnSH{ AM8*( R-E) 1 
CHBPE=COSH( A!^B*( R+E.n 
CSBHE=COS(AYD*( B-E) } 

SN8PE=SIN{ A*1B=^(B+E) ) 

SHBME=SINH( AHB*{ 8-E) ) 

SHBPE=SINH( AP'B={'C B+E) ) 

SNBME=SIN{ AMB*( B-F) ) 

EM=1.0/ (PI*SQRT< 2.0) *THETA*{ SHE^SHE- S‘4E’»‘SME ) ) 
FH1=2.0#.SHBPY*SN8PY 

TW1=EH1*{ SHE SR PE *CHBME-SNE *CHBPE=i=CSBH E ) 

TW2=CHBPY*SNRPY~SHBPY*CSBPY 

TW3=SHE*(SNBPE*CHBME-CSBPE«SHBME) 

TW3»TW3 + SNE* CSHRPE^CSBME-CHBPE^SMBME ) 
TW=TW2*TW3 

AEW(I,i)=EM’^={TWl+TW) 

U=0. 15 

0P=1.0+U 

0M=1.0-U 


TP=3.0+U 

PHI=PI*Y/B 

XI=Pl-ABS{PHI-Sn 

TX=THETA*XI 

SHX=SIMH(TX) 

CHX=COSH{TX) 

T$I=THETA*SI 
THI = THETA*PHI 
CHI=COSH(THI ) 

SHf = SINH(THl ) ■' 

C HSI »COSH { TS T ) - ' ■ ' "’r ' ■' ■ My/.; 

SHSI=SINH{TST) ' 

EM=-l.f)/f4.0*SIG*SHG’)‘SHGI y' 

TWl= (OM#SIG-TP*SHG) THl ♦SHG^'SHI *0M J •I'f 

.'.TWlsTWl/ {TP*SHG*CHG-0»A'SfSIG> 


T Wl= TWl * ( ( OM’S'S .IG^CHG-OP^SHG 1 *CHS I -3'^4Ti|*SWG« S»SI' i''U' "iff 
TW2=*STG#CHG^SHI=«>0M-0W*THI=»'SHG*CHI '■ 

TW2=f TW2# ( ( 2» 0*SHG+S| G*CHG*OM J *SHS'l^-T'SI*«^SHG#tHS I *0H 

tw 2 =tw 2 /(tp*shg*chg+o«#sig)i MM/'-v;::"'’-'. 


' p .jfl!,.''’ ^ 


V iv; 
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TW=TW1+TW? 

TW= TW+ ^OMsi'S JG-i^CHG-TP^SHG) =i'CHX-rx«=SHG*SHX=:'.T''' 

AEU(I,J) = Ef'1*TW 
10 CONI’' '‘JUE 
DO 15 T=ls3 
DO 15 J=l,3 

WtEW(;i,j)=Al=y (T^ j) 

15 MEU{I,J)=AeU{T, J) 

DO 20 J=l,3 
MEU( J,4) = ARU C J, 2) 

20 WEWt Jf4) = AEW ( J,2 ) 

RETURN 

END 

SUBROUTINE GFN{A,AA> 

# # 5^ ^ 3^,: j|c ^ ^ ric ^ «|s t -K -Y 


THIS EXPANDS THE KOjAND KJ MATRICIES OF 5X9 TO 9X^ j*:! .JG 

MAXWELL S RECIPROCAL THEOREM 


DIMENSION A{5,9),AA19»9) 
on 10 1=1,5 
DO 10 J=l,9 
10 AA{I,JI=ACI, J) 

DO 20 1=6,9 
DO 20 J=l,9 
IF(J.GE.6) 60 TO 30 

AA{I,J)=AA{K,n 
GO TO 20 
30 L=10-I 


K = J-4 

• AA(I , J)=AA(K,L) 

20 CONTINUE 
RETURN 
END 

SUBROUTINE DIST RI ( A, SUL, KO, KUI 


-i*. <?' 




THIS GENERATES TWO 5X9 MATRICI'ES S' ONE 'EAC'H FOR XO AND' K i ''"'g 

, ■ ‘:r.» 

, ■ .... ' .. ;■ '' 

DIMENSION K0{5,9),KUf 5,9) . ' ' ^ ,,,v,. -v ■ 

REAL KO,KU,MEW,MEU 
THETA=A " 

PI=4.0#ATAN{ 1.0 : : 
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B = 1*0 
y=-1.25*B 
SI G=THETA’!‘PI 
CHG=COSH^ STG) 

SHG=SINHCSIG) 

DO 10 1=1,9 

E=-0.25*B 

r=Y+0«25*B 

AM8=PT^THFTA/{ B*S0ST{2.0) ) 

TME=2.D*AMB*R 

SHE=SINH?TyE) 

SNE=STN(TM=) 

DO 10 J=l,5 
E=E+0,25*B 


IF(I.EQ.6,AND.J«EQ«l) 

TFCI,E0.7. AMDe J,LT„3) 
TFCI.FQ,P,AMn.J.LT*4) 
T F(I .EQ,9.AND.J.LT.;.) 


GO TO 100 
GO TO 100 
GO TO 100 
GO TO 100 


30 CONTINUE 

CH8PY=C0SH{ AMB*C e+Y) ) 

CSBPY=COSt AMB*(B+Y) ) 

SNBPY=SIN{ AMB*(B+Y) y 
5HBPY=SINH{ AKB’i'C B+Y) ) 

CSBPE=COSCAMB=^fB+P)) 

CHBME=COSH{ AMB#{ B-E) ) 

CHBPE=COSH{ B+En 

CSBME=COS( AMB*{B-E> ) 

SN8PE=S IN^ AFie*(B+F) ) 

SH8ME=SINH{ B-E) ) 

SHBPE=SINH( AWB*{ R+E) ) 

SNBME=SIN(ANB*(8-E) ) 
E«=THE/(SHE=«'SHE-SME'=^‘SNE) 

EM1 = 2.0*CHBPY’S'CSBPY 

TW1=EM1*CSHE*CSBPE*CHBM E-SNE*CHRPE=8'C SBW E > 

TW2=CHBPY*SNBPY+SHBPY*CSBPY 

TW3=SHE=>CSN8PE*CHBMe-CSBPE’5'SH8ME) 

TW3=TW3+SHE*{SHBPE’«‘CSfiME-CHBPE*SMBHE ) 

TW=TW2*TW3 

K01U,I)=EM*{TW1+TW) 

PHI=PI*Y/B 
SI=Pl=«'E/8 
XI=PI~ABS{PHI'-SI ) 

TX*THETA=i'XI 

SHX=SINH?TX) 

CHX=C0SHITX1 
TSI^THETA«SI 
TH1=THETA=^PHI 
CHl = COSHfTHf ) 

SHI=S1NHITH1) 

CHSIs=COSH(TSn , 

. f- ■" 



oooooooo 
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SHSI = SINH (TSI ) 

'?PHI = f S!G^CHG-SHG) ^CHI-THI*SHI*SHG 
^SI={S!G*CHG“SHG )4:CHS|-TST*SHSI*SHG 
0PHT = t2,0*'^HG+SIG^CHG)^SHl“THl=i'SH0>«'CHI 
0S!={2,0#SHG+SI G=«‘CHG)^SHSI-TSI^SHG>«'CHS J 
rM=SIG/{2«0*SHG*SHG} 

KU( J,I)=FM=^'{ f SJG>J-CHG+SHG) *CHX-TX *S HG^'S HX +R P HT / f " , . -J' : i-G’?- C MG- ' 
1) +QPHT*QSI/( ?.aX'SHG«‘CHG+S IGH 
T Fd.EQ.e. AND* J«E0« 1) GO TO 110 
TFa.F0.7»AND* J*LT. ?) GO TO 110 

TFd.EQ,R»ANn« J*LT.i) GO TO 110 

TF(I.E0*9.AND. J.LT.:i) GO TO 110 
GO TO 10 
100 YI=Y 
EI=E 
^Y=-Y 
E=-E 

GO TO 30 
110 E=EI 
Y = YI 

10 CONTINUE 


RETURN 

END 

SUBRCUTINE UNCON (N, X, Gf E, IMAX»R) 


4 : » * * * 


♦ 4 ^ ♦ ^ 


THIS INCORPORATES THE D.F.P. ALGOPITHY FOR UNCONSTRAINED HINIYI- 
-ZATION OF PENALTY FUNCTION F 

C0FM0N/SETl/SPAN,WTFjWLFfNBiTH£TA,T!>L,PL,»,WP70R,WPFP|?, B«FLL»B*^Frj 
lL,BMWC»BMLL,FPTWP,FPTTP,FPBTPt FPBWPr SCO, FS , CC, CHTW, CAC, GKS, RE 
2,C1TW,CS 
C0NM0N/SET2/ITER 
C0NM0N/SET3/ PT, POfNTCjNC, APC 

DIMENSION X{10),Ga5J» GR ADUO ), ST 10 ) , GADP f 1 0) » AH( lOt 1 0 1 r 

1YQ(10),HY(10) 

ZjXDClO) 

dimension TPL{4) ,PL{7)»PC71,WP70R,C9)»»<PFP{9) 

5 CALL FTN{G,X,F,OBJf M,N, R) , ' 

ITER=1 " .. ■ ' 

T0=0.10 
DO 10 I»lf M 

IFCG(l).GT.O.O) GO TO 160 ; 

10 CONTINUE 

15 CALL GRAON(N*X,MyR,GRADl . . 

20 DO '30 I«lf N 

DO 30 J=l,N ■ . .. . '-v ■ 
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&HCI,J)=C«D 

TFH.EO.J) AH(I,J) = 1,0 

30 CONTIMUE 

40 suM=a. 

DO 60 T=1,N 
A = 0, 0 

DO 50 J-ltN 
50 A=A~AHCI, JI^GRADCJ) 

S!JM=SUM+A*A 
60 St I )=A 

SUM=SORT( SUM) 

DO 65 I=1»N 
65 StI)=S(!)/SUM 
I K=0 
A”0#0 

no 70 1 = 1 , N 
GADPtI)=GRAO{I ) 

70 A=A+S(T )*GRAO{I ! 

IFtA.GT.O.O) GO TO ?0 

IFtABSCA)«l.T»0,5S-035 GOTO 150 
no 15 1=1, N 
75 XDtI)=X{I) 

F1=F 

CALL FI8B0(X,N,S,F,iMAX,T0, AL, R,M) 

79 CALL FTN(G,X,F,08J,M,N,R) 

CALL GRADN(N,X,M,R,GRAD) 

PRINT 601, ITER 

601 F0RMATt5X,«ITER =^,I3) 

PRINT 602,{XtI),T=l,H) 

602 FORMATtSXf^VARIABLFS ARE*, / , 8E 15 . 7 ) 

PRINT 603,{G(I),T=1,F) , 

603 FORMAT{5X,*CONSTRAINT ARE*, /, 9E12, 4, /, 6 El 2« 4) 

PRINT 604, F, OBJ 

604 F0RMATC5X,*F =*, E15,8,*0BJ=*,ei5,8 ) 

PRINT 5i5,PT,P0,NTC,NC, APC 

515 F0RMAT{2E20,8,2llf>,E20.8> 

1TER=ITER+1 ^ ■ 

IFtlTER.GT.lO) GO TO 150 

77 EU=0.0 
El=0. 

DO 78 1 = 1, M '■ ■ 

A=XD(I)~X(I ) , , ' . ^ 

|F(A.GT*eU) EU=A , vvP"v.;''' '' 

E! = EI+GRAD{n*GRAO(n *''■■■' ; ' '' i 

XD(I)=X{I) 

78 CONTINUE 
,EI=SQRT{'EI) 

F0?=(F1-F)/F 

IF(ABS(EU),LT. 0.005. OR. ei,LT.O*OOOl.OR.,EO,L;T.O*OOOII TCl,,' 150 

ri -»C '-4. /i , > , . , ‘ ' ’ ' ' T • y- i, > 'l". , s'-f ,,, ' > /t, 
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do 80 1=1, M 

80 YOa}=GRAD(T)-GADP(I) 

STY=0,0 
DO 90 !=1,N 
90 $TY=STY+S{n’!'YQ{I) 

TF(STY«L!:.0, 1E~155 00 10 40 

DO 110 1 = 1, N 
A = 0.0 

DO 100 J=1,N 
100 A = A+AHnf J)*Y0( J) 

110 HY(I)=A 
YHY=0.0 

no 120 1 = 1 , N 

120 YHY=YHY+HY(| )*Y0n) 

DO 130 1 = 1, N 
DO 130 J=1,N 
AM=StT)=«‘S( Jl/STY^S'AL 
ON=-HY(I)*HYCJ)/YHY 
130 AH(I,J)=AH(T,J)+aM+ON 
GO TO 40 


150 RETURN 

160 PRINT 170, CXCI ) , I=L,N) 

170 F0RMAT(5X,=^STARTTMG POINT VI OL AT ES DONSTRAINTS >«'/ , 8E I 5 • 6 ) 

PRINT 180, {G(I), 1=1, N) 

180 F0RMAT(9E12»4,/,4F12»4) 

SO TO 150 
END 

SUBROUTINE GRAON( N, X, M, R, GRAD) 

THIS GRADIENT OF THE FUNCTION AT A POINT K USING CENTRAL DIFFEReNC 
SCHEME 


C0NM0N/SeTl/SPAN,MTF,HLF,NB,THeTA,TPL,PL,P,WP70R,WPFP,B,BMFiL,B^FD' 
1L,BMWC,BMLL,FPTWP,FPTTP,FPBTP, FPBWP, SCU, FS,CC,CHTW,CAG,SK$,RE ' 
2,C1TW,CS 

dimension TPL{4) ,PL(7},P(7),WP70R(9),WPFP{9) 
dimension XdO) ,G(lS),GRADf 10),XP< 13 ),FPUO},FH{ 10) ,, 

DO 5 1=1, N 

5xp(i)=x{|} ■ 'Y"V' ; 

0010I = 1,N '■ <«- 

xp a )=x(n*i«oi . 

CALL FTM{G,XP,FF,OBJ,M,N,R) . , '' 

xpd )=X{I ) 

10FP{I)=FF ■'/ -■ 'rr'’"*" '*■' 'f ' " ,>■ ■. ■ :■ ■ 

DO 20 1 = 1, N ^ 

XP{I) = X{I )^0»99- ' '■ ' ,1 r • - , 
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CALL FTNfCt XP? F F, OBJ, N, R ) 

xp(i)=xn) 

20 FMd )=FF 
00 30 1=1 ,N 

30 SRAD{I)=(FP{I)-FM(I} )/{0,02*X{ I) ) 

RETURN 

END 

SUBRCUTIME FT BBD(X, M, S, £, TDt AL, R, M ) 

:l)c ^ sie ;;jc ^ 4; # 5{c ^ jjc ^ jji: sj; ){j ^ 5^ 5>- ^ ^ j- 


THIS FINDS THE OPTIMAL STEP LENGTH IN DIRFCTJON S BY iifX' 

METHOD 


3^ ^ ^ ^ ^ ^ ^ ^ 5§i ^ ^ 3|S 5^ 5^ ^ ^ ^ 5^ J.. 

C OMMON/ SETl / SP AN, W^F sI^LFjNB, THETA, TPL,PL,P, WPT OR,^^PFP, P, F’MFLL, P'^FD 
1L,BMWC,BMLL, FPTWP, FPTTP, FPBTP, FPB4P, SCU, FS, CCj CHTWtCAC s GKS, RE 
2<C1TW,CS 
COMMGN/SETP/ITER 

dimension X(10) ,G(15), FI( 20) ,TPL( A), PL{T), P( 75 jrtPIO*’ { 91 t WOFP(H) , 
IXAdO) ,XB(10 ),S( 10) 

AL=0.0 

KJ=0 

STEP=0»0 

XMAX=15 

K=0 

TF(R.GT.5«0) TO=U()/FLOAT( ITER) 

IF(R.LT«5«0) T0=0,5/FL0AT( ITER) 

DO 5 1=1, N 
5 XA(I)=X{I ) 

CALL FTN(G,XA,FA,OBJ, M, N, R) 

10 DO 15 1=1, N 
15 XB(I)«XA(I)+Tn*S( J ) 

CALL FTN(G,XB,F8,0BJ,M,N, R) 

IF(FB,LE.FA) GO TQ 20 
GO TO AO 
20 DO 25 1=1, N 
25 XA{I)=XB(I) 

STEP-STEP+TO 

FA=FB 

K=K+1 . , .. 

SO TO 10 ’ 

AO A 1=0,0 , / . ' 

B1=TC ■' 

, DO Al 1 = 1, N - ■■ , ^ '■ ■■■:;' ’-p : 

4l’'x(I)=XAlI ) ■ H 

A2 FI(l)=1.0 ' 

FH2)=1,0 

' DO AS 1=3, 1 max, , ^ 
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45 F! a) = FI {I-2) + FMT-l) 

!!=!MAX 

AL2S = Fn JT“2)/Fnn 
J = 2 

900 AL1=B1-A1 

!F{AL?S.GT,AL1/2.0) GO TO 300 

T FUL7S,ee, ALl/ 2 , 0 ) GO TO 111 

X1=AL2S+A1 

X2=B1-AL2S 

GO TO 400 

111 Xl=Bl~AL?S--0.005 
X2=Al+AL2S+0o005 
GO TO 400 
300 X1=B1-AL2S 
X2=AL2S+A1 

400 DO 401 I=i,N 

XB(I )=X(1 )+Xl*S( I) 

401 XA{| 1 = X(T ) + X2=^S( I) 

CALL FTN{G,XAfFlsOBJtM,N, R) 

CALL FTN{G5XP,F2tOBJ,M,N, R) 
1[F(F1-F2) 600f4‘"0,500 

450 A1=0,0 
B1=X1 

IF(J.EO« (N-l }) GO TO 700 

N1=N-J-1 

N2=N-J 

AL2S=FnNl)/FI (NZI^iBl-Al ) 

J = J+1 
GO TO 700 
500 A1=X1 

AL2S=X2-A1 
GO TO 700 
600 B1=X2 


AL2S=X1-A1 
700 J*=J+1 

IFCJ.EO.N) GO TO 800 
GO TO 900 
800 STEP=STEP+A1 
AL=AL+STEP 
RETURN 
END 

SUBROUTINE FTNtG,X,F, OBJ, M,N,R3 ' 


THIS FINDS THE PENALTY FUNCTION, OBJECTlYf, 'FUMCTfON AMO CONSTCAl'l 
TS' FOR BOX TYPE SECTION AT A‘'0ESIGN POINT X-f. ^ '-j-S' ' ' 

subroutines REOUIREO' ARE,^ MAlN,GEM,O|ST;:RliTIJAG0EfUNCfmirSPADM''^^ 

FIBBO . '■ ^ . V V ’4 



c 

•£ ^ ^ ^4: s^c ^c ?»j: it ^!s ? : ;}' ;■': 

C0MM0N/SFT1/£PAN,HTF, WLF,NE, THETA, TP L, PL, p, WP" 0^ , F P , F , •■ 'TL L, ^''F' 

' 1 L jBMWC, BMLL, FPTWP, FPTTP, FPBT P, FPBWP, SCU, FS, CC, CHTw, C^C , GF S L" 

2,CITW,CS 

COMMGM/SET3/PT,PO,NTC,MC, APC 

niMENSlON X( 10),GC1!: ),KQ(5,9),KJ( G,9),MEU( 3f , A ), 'O , 

IWPFP (9) ,0K(9,9) , ABKO( 9) , ABKU{9),TPL{ A), CU{ 6 ), PL ( ■’) , Pf ~ 1 
PEAL KOjKUf PFH, MEU 
MTB=25.0 

-IU*XCA)-X(1)-X(2) 

APC=WLF*(X{i )+X(2) )+HU»2.0=«X{3) 

FMPC8=WLF«( X(2)AA2)*0.5+HU’«!X?3 )*{X( 2 )+HU=»=:*f )*?« "'+WlF*/{ j )=F{ xf A) ~ 
1.55^X{13 ) 

YPB=FHPCB/APC 

yPT=Xl43-YPB 

C^J = WLF=«‘(X( 1 )^*3 3/l?®a+WLFAX{l C YPT-D,5*X( U )^^?) + 2.. ' * 4(3) 

13) /12«r)+2,0«X(3) { YPT"HU*0.5-XC 1 ) ) **?) +X ( 2 ) F /I 2. 0+ULF4«Vj 

22 )*( (YPB-0,5*X(2) )**2) 

ZPB=CMI /YPB 
ZPT=CMI/YPT 

HA=HU*(WLF-2»0’i'X (3) ) ' 

Dl=4.0*HA=f^HA/( (WLr~2.0*X{ 3) )/X( 1) + {WLF-2.0*X(3) )/X( 2) + 2«0><'HU/X(''») ) 
T J=WTF=«' ( X ( A ) ) /1 2 . 0“ ( WT F-2 . 0=(‘X ( 3 ) ) HU**3 ) / 1 ?« 0 

HT={ WTF-?.0=f'X{3) )=«'HU , 

0J=4,0*HT/{ (WTF“260=<^X{ 3) ) /X( l)+{ WTF- 2, 0*X { 3 ) ) / X { 2 ) +2« 0«:HtT/X C 3) ) 
TMB=WTF*(X(2)A*2)=4'0.5 + HU*X(3)=S=(X( 2) + HU>!'0.f )=*<2*0+HTF*>:X( 1 ) *«•( XCA ) -0,5 
1*XC1)1 

ATB=WTF*CX(1)+X(2) )+HU=5'2.0=«'X{3) 

0T8=TM8/ATE 

DTT5:X(A)-[^TB 

ZTBB=TJ/DTB 

Z.TBT=TJ/OTT 

SI=CNf/WLF 

SI0®0!/WLF 

SJ=TJ/WTF 

SJO=OJ/WTF 

THETA*8*{ (SI /SJ) **0,25) /SPAN 
ALPHA=0.2175*(SIO+SJO)/SQRt{ SI*SJ) 

SALPH=SQRT(ALPHA) 

CALL DISTRICTHET A,ALPHA,KO,XU) 

CALL GENCKO, OK) . 

DO 10 1=1,9 
XX=0,0 
00 5 J=l,9 

5 XX=XX+WPTOR( J)*OKII, J) 

10 ABKO(T)=XX/A.O 
CALL GENCKUfOK} 

DO 20 1=1,9 

XX=0.0 - V 
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DO 15 J=l,9 

15 XXsXX+WP^'ORC J)*nK{I,J) 

20 'XBKUd )=XX/4.0 

Y=0.0 

no 25 1=1,9 

XX=ABKU(I )-ABKQd) 

ABKO(T) = 4PK.O(n+XX*£ALPH 
IF(Y.,GT. APKC (I) ) GO TO 25 
Y=ABKO<n 
T.J=I 

25 CONTINUE 
TF>=IJ-1 

nlFB=FL0AT(I^'>*B/4«0 

XX=WLF/2,0 

MS=N8-1 

DO 30 1=1, NS 

MBS=I 

Y=DIFB-XX 

IF(Y. LT.WLF) GO TO 40 
30 XX=XX+WLF 

40 0KW=ABK0(IM)+4,0>!'Y4( ABKOC !J)-AB<D( I'^n/B 
CALL GEN(KO,OK) 

DO 45 1=1 ,9 

XX=0.0 

DO 42 J=l,9 

42 XX=XX+WPFP{ J J) 

45 ABK0(I)=XX/2,0 
CALL GEN (KUfOK) 

DO 50 1=1,9 
XX=0.0 
DO 48 J=l,9 

48 XX=XX+WPFP{ J) 

50 ABKUCI)=XX/2«0 
DO 55 1=1,9 
XX=ABKU(1 )-ABKO( n 
55 ABKO<n=ABKO(I) + XX=«'SALPH 

DKF=ABKO{IM1+4.04{ ABKO( U )-ABXOC IM) J/B 
DHD= C WLF’I'X (4 )-APC)*WTB*0. 0024 
WS=APC*0.0D24 

B MS=WS4SPAN*SP AN/8, O+DWD’S'SPAN/2,0 
BN=FLnAT(NB) 

PBMFLL=1.1*DKF=«^BMFLL/BM 
PB!«FDL=l,l*DKF’^BMFDL/BN 
PBMWC=BMWC/BN . 

PBMiL=BMLL*1.141,15»!'DKW/BH :■ /; ; 

PBMA=PBMFLL+PBMFDL+PB«WC+PBMLL ::T/G::4'f' 

PBMD=8MS+PBNFDL+PBF!K’C - 

PBMi=PBMFLL + PBMLL . ‘''T ■ , 

2TP= (PBMA + ( 1.0-1. 0/Re)*BMSi/1FPTMP-FP'TTP/Rf'l,.^ 
’ Z 8P= ( MA+ ( 1 .0-1 • 0/RE 1 *BWS ) / ( PPBTP/R WPBWP P 'T? 



3(1)=ZTP/ZPT-1.D 
G ( 2) = ZRP/ZPB'-l.C.t 
pp=:YPP-'),08R=^X{4) 

i>0=( PRHA + P'‘'tS + FPBWP=f‘ZPe)’i'RE/( EP + ZPB/APC} 
FPBTA=Pn/APG+P0*EP/ZP8-BMS/ZPR 
G(3)=FPBTA/FPRTP-1,0 
FPTTA=P0/APC-P0*FP/ zpt+BMS/ZPT 
.G(4)=FPTTP/F PTTA-1, 0 

FPBWA=(Pn/APC + PO*FP/ZPB)/RE -BMS/Z PP-PRF. A/ZPB 

GC5)=FPBWP/FPRWA”l#a 

TF(ABS<G{5) )<.LT.O,Of 1) G(5) =0.D 

FpTWAs{Pn/APC-PO*FP/ZPB) /RF +BMS/ Z PT +P BM A/ Z PT 

G (6) =FPTWA/FPTWP-I .3 

CNi=PC/{0.6*FS*0«38 5>i=12.0) 

XX=AM0D(CNrl«0) 

IFCXX.GT.O.O ) GO TO 60 
MC=IF!XCCM) 

Y=0.0 . 

AOC=0,0 
GO TO 65 
60 Y=1.0-XX 
ADC=12.0*XX 
CN=CN+Y 
NC=IFIXCCM) 

65 ATS= {(CN-Y)^12.D+ADC)*0.385 
fFD=YPT+FP 

UMR=1.5*PBMA+2.5*PBML 
ATSF=0.6P=«‘SCU^( WLF-XC2) )*X{1)/FS 
ATSW=ATS~ATSF 

G(10)=ATSW*FS/CWLF*E?:D*SCU*0*24)-1.D 
USM= ATSW*FS=^EF0*(l«0”0»75*ATSW*FS/f X { 2 ) ’«'eFDYSCU) ) 
IFCXll) .GE.EFO/4.0) GO TO 70 
USW=USM+0,7’«'SCU*X( 1 )*(hLF~X( 2) )’«=( EFD-0,5*X{ 1)5 
70 G(7)=UMR/USf^-1.0 

CALL TRACOEfTPL, MEU,MFW, THETA) 

00801 *lf6 

XX=0.0 

!F(I,GF.4) GO TO 72 
DO 59 J=l,4 
59 XX=XX+MEWI If J5 
EU!I) = XX 
GO TO 80 
72 XX=0.0 

T0=I-3 ' 

DO 77 K=lt4 
77 XX=XX+MEU{I0,K5 

Eucn=xx 

80 CONTINUE 

DO 82 I = lt5,2 , 

E u ( I ) »E u ( 1 5 + 1 eu c 1 + 1 5 -Eu { I ) I *s At PH. . : 



82 EU(I + l)=n,0 

P I =4, 0* AT AN ( 1,0) 

YML=0,0 

TK=0 

DO 90 1=1, 5, 2 
f K=I K+1 
FI=FLOAT{ J) 

XX=0.0 
DO 83 J=l,7 

85 XX=XX+P( J)*S?N{EI*P!*PL( J)/SPAN) 

IJ=IK+1 

YML=EIJ{ I)*XX^( 

90 CONTINUE 

YML=YML*?s04B/SPAN 

TYMl=YNL*UTF 

PT=(TYML+FPBWP=«'ZTRB+FPTTP*ZTBT/RE>TATB4RE/{ 7TBB+ZTRT ) 
G{8)=PT/ATB-1.0 ■ 

S{ 9)= {ZTBT4PT + Re*TYML)/{RE*ATB*ZTBT=!<FPTWP)-i,0 
PT=5.0=<'PT 

TNC=PT/(0.6=^FS*0.?.83^12.0 j ’ 

XX=AKODtTNC, 1*0) 

■ IF(XX,EQ,0«0) GO TO iOO 

Y=1,0-XX 
TNC=TNC+Y 
100 NTC=IFIX(TNC) 

3(11)=15.0*X(4>/SPAN-1«0 
GC 121=14.0/X(1)-1.0 
G{ 13 ) = 14.n/X{2)-1.0 
G{14) = 14.0/X(3)~l«0 

XX=( APC=«'SPAN + (X(4)-X{ 1 )-X( 2) ) *( WLF-X C 3 )*2* 0 ) 4^. 0*WT8 ) «CC 
XX=XX+C FLOAT(NC) *CS + ( ATS’S'ClTW/12.0) )*SpAN4{ ]. + 8 . *F p*F; P /{".<« SPAN^'SP 
IAN)) .. '■ 

0BJ=XX*BN+TNC*CHTW*2,0=«'8 . . 

D B J= OB J +T NC* C A C + CN*C AC* BN 

SKS=0.0 

DO 110 1=1, F 

IFCG(I).F0.0«0) GO TO 110 
GKS=GKS+1.0/G(I ) 

110 CONTINUE 

F=OBJ-R*GKS 
DO 119 1=1, M 

IF(G(I) .GT.O.O) F=0.5E 11 

119 CONTINUE 
RETURN 
END 



